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ABSTRACT 


Several  autonomous  satellite  navigation  configurations  were  studied  on 
a  common  low-altitude  basis  with  the  emphasis  on  accuracy  and  convergence. 
Analysis  was  performed  by  simulations  using  recursive  filtering  techniques. 
Configurations  investigated  included  horizon  sensors,  unknown  landmark 
tracking,  known  landmark  tracking,  and  one-way  doppler  systems.  Relevant 
sensor  combinations  were  analyzed  to  determine  possible  performance 
augmentation. 
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1.  INTRODUCTION 


This  report  documents  a  four-month  study  on  autonomous  satellite  navigation 
in  which  various  sensors  were  compared  on  a  uniform  basis.  Special  atten¬ 
tion  was  given  to  possible  hybrid  configurations  for  optimizing  autonomy, 
accuracy,  simplicity,  and  time  of  convergence.  A  special  computer  program 
called  SNAP  (Space  Navigation  Analysis  Program)  was  developed  for  analysis 
of  the  various  systems.  Comparisons  were  made  on  the  basis  of  simulation 
results  for  a  low  altitude  satellite  orbit.  Among  the  sensor  types  considered 
were  unknown  landmark  tracking,  known  landmark  tracking,  horizon  sensors 
and  one-way  doppler  systems.  The  work  of  previous  investigators  in  these 
areas  (Refs,  i,  5,  6,  7,  20,  21,  22,  23,  24,  27,  33,  35,  and  41)  was  used 
in  developing  appropriate  system  models  and  in  choosing  the  parameters  for 
the  cases  considered. 

The  conclusions  obtained  in  this  study  are,  of  course,  influenced  by  the  specific 
configurations  analyzed  and  the  low  altitude  orbit  used.  Hence,  care  should  be 
exercised  in  extrapolating  the  results  reported  herein  to  other  situations. 
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2.  SIMULATION  SOFTWARE 


The  various  sensor  configurations  considered  in  this  study  were  evaluated  by 
means  of  processing  simulated  sensor  data  via  a  suitable  estimation  algorithm. 
The  Aerospace  Corporation  has  extensive  software  capabilities  for  performing 
such  analyses.  Among  the  available  programs  are  TRACE  (Refs.  9  and  42), 

MVS  (Ref.  39),  ONUMS  (Ref.  1)  and  ONAP  (Refs.  25,  26,  and  36).  The 
latest  program,  SNAP,  was  specifically  designed  for  this  study  to  facilitate 
simulation  of  hybrid  sensor  systems.  The  basic  simulation  logic  used  is 
shown  schematically  in  Fig.  2-1.  Some  of  the  principle  features  of  this  pro¬ 
gram  are  highlighted  below. 

2.1  By  means  of  a  simple  input  flag,  the  user  can  choose  between  Standard 
Kalman,  Stabilized  Kalman  or  Andrews  square -root  filtering  algorithms 
(Ref.  2  and  18).  The  Andrews  filter  is  a  generalization  of  an  approach  orig¬ 
inally  suggested  by  Potter  (Ref.  3).  Dynamic  "forcing  noise"  can  be  accom¬ 
modated  with  square-root  filtering  by  use  of  either  a  Cholesky  or  a  House¬ 
holder  decomposition. 

2.  2  In  order  to  simplify  programming,  all  partial  derivatives  are  computed 
numerically.  While  analytic  partials  are  generally  more  desirable  in  a  specific 
application,  the  numerical  technique  was  used  to  provide  flexibility  as  simu¬ 
lation  demands  changed.  In  a  limited  number  of  test  cases,  the  results  obtained 
using  this  approach  compared  quite  well  with  the  results  of  other  software  which 
was  based  on  the  analytic  approach  to  partial  derivative  computations. 

2.  3  The  dynamic  models  to  be  used  for  any  run  can  be  chosen  from  any  of  the 
following  "real  world"  or  "filter  world"  options: 

Real  World  Models  (See  Appendix  A) 

Keplerian,  no  atmosphere 

Keplerian,  U.S.  Standard  Atmosphere,  1962 

APL  "8-8"  geopotential,  no  atmosphere 

APL  "8-8"  geopotential,  U.S.  Standard  Atmosphere,  1962 
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Figure  2-1.  Simulation  Procedures 
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Filter  World  Models  (See  Appendix  B) 

Keplerian,  no  atmosphere 
Keplerian,  exponential  atmosphere 
p  and  geopotential,  no  atmosphere 

p  and  geopotential,  exponential  atmosphere 

These  models  were  chosen  primarily  to  provide  a  reasonable  mismatch 
between  the  real  world  and  filter  world  emphemeris  (Refs.  5,  6,  21,  22, 

23,  24).  Extensive  investigation  would  clearly  be  required  to  determine 
optimum  filter  world  models  for  an  operational  system. 

2.4  A  simple  array  of  integers  specified  at  execution  time  is  used  to  establish 
the  elements  of  the  state  vector  to  be  estimated  for  a  particular  case. 

2.  5  Configurations  for  autonomous  navigation  can  be  chosen  from  among  the 
sensor  options  listed  below  by  an  appropriate  array  of  integers.  Up  to  three 
independent  sensors  of  any  specific  type  can  be  used  simultaneously,  with  the 
number  of  different  types  chosen  for  any  configuration  limited  only  by  the 
rationality  of  the  final  system.  Furthermore,  an  independent  measurement 
schedule  can  be  provided  for  each  sensor  chosen.  As  required,  any  sensor 
observation  can  be  biased  with  or  without  that  quantity  appearing  in  the  state 
vector. 

Sensor  Options  in  SNAP 

1.  Azimuth  and  elevation  to  unknown  landmark 

2.  Horizon  sensor  (direction  to  center  of  earth,  angle 
subtended  by  disk  of  earth) 

3.  Azimuth  and  elevation  to  known  landmark 

4.  Range,  azimuth,  and  elevation  to  known  landmark 

5.  Range  to  known  landmark 

6.  Range  rate  to  known  landmark 

7.  One-way  doppler  to  known  landmark 

8.  Star  elevation  above  earth  horizon 
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9.  Range  to  satellite 

10.  Azimuth  and  elevation  to  satellite 

11.  Range,  azimuth,  and  elevation  to  satellite 

12.  Range  rate  to  satellite 

13.  Range  difference  between  two  satellites 

14.  Range  rate  difference  between  two  satellites 

15.  Range  difference  between  two  known  landmarks 

16.  Range  rate  difference  between  two  known  landmarks 

17.  Radar  altimeter 

Since  the  total  time  available  for  this  study  was  quite  limited,  no  attempt 
was  made  to  evaluate  all  of  the  above  sensor  systems.  Specifically,  only 
options  1  through  7  were  examined  in  this  study. 
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3.  BASIC  ASSUMPTIONS  AND  PROCEDURES 


In  order  to  provide  a  uniform  basis  of  comparison  for  the  systems  considered, 
the  following  basic  assumptions  and  procedures  were  used  in  all  the  simulation 
studies: 

3.  1  The  basic  simulation  procedure  used  for  all  studies  is  shown  schematically 
in  Fig.  2- i  . 

3.2  A  reference  orbit  having  a  nominal  80  n  mi  perigee  and  a  200  n  mi  apogee 
was  used.  This  orbit  was  arbitrarily  chosen  as  typical  of  low -altitude  orbits. 
Detailed  orbital  conditions  are  given  in  Appendix  C. 

3.  3  Unless  otherwise  specified,  real  world  data  were  generated  using  the  APL 
"8-8"  geopotential  and  the  1962  U.S.  Standard  Atmosphere  models  described  in 
Appendix  A.  All  filter  calculations  were  based  on  the  p  -  geopotential  and 
the  exponential  atmosphere  models  described  in  Appendix  B. 

3.4  Since  the  total  time  available  for  the  study  was  quite  limited,  no  attempt 
was  made  to  compensate  for  dynamic  model  errors.  Where  possible,  however, 
an  approximate  point  of  divergence  was  identified.  Note  that  the  basic  problem 
of  compensating  for  model  mismatches  has  been  shown  to  be  solvable  empiri¬ 
cally  (Refs.  6,  22),  generally  requiring  extensive  analysis  and  simulation 

for  a  specific  configuration. 

3.5  Because  of  recent  advances  in  attitude  reference  systems  (Refs.  10  and 
11),  attitude  errors  were  assumed  to  be  negligible  compared  to  navigation 
sensor  errors. 

3.6  A  spherical  earth  was  used  for  all  geometric  calculations. 

3.7  Unless  otherwise  specified,  the  Andrews  square-root  filter  was  used  for 
data  processing  to  reduce  possible  roundoff  errors. 
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3.8  A  specific  format  for  presenting  the  results  of  each  system  studied  was 
used.  First,  a  detailed  description  of  a  baseline  configuration  is  presented  in 
tabular  form.  Then  results  of  simulating  the  given  system  are  discussed. 
Deviations  from  the  baseline  configuration  are  explicitly  noted  and  the  corres¬ 
ponding  simulation  results  are  discussed  in  turn.  The  systems  are  evaluated 
in  terms  of  a  root-sum-squared  (RSS)  position  error;  graphical  time  histories 
of  position  errors  are  given  for  several  key  cases.  Most  of  the  data  obtained, 
however,  are  summarized  in  tabular  form. 
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4.  UNKNOWN  LANDMARK  TRACKING 


With  an  unknown  landmark  tracking  system,  satellite  navigation  is  accomplished 
by  making  repeated  measurements  of  the  inertial  direction  from  the  satellite  to 
a  point  of  unknown  location  on  the  earth's  surface.  This  so-called  "unknown 
landmark,  "  is  visible  to  the  satellite  for  only  a  relatively  short  time.  After 
visibility  is  lost,  another  suitable  unknown  landmark  is  selected  to  provide 
additional  navigation  information. 

This  type  of  system  has  been  extensively  investigated  in  the  past  (Refs.  4,  5, 

12,  14,  21,  22,  23,  and  32).  The  most  recent  study,  by  Hendrickson,  shows 
that,  in  spite  of  certain  numerical  problems,  long  term  stable  operation  can 
be  achieved  if  the  data  is  processed  with  either  a  square  root  algorithm  or 
with  the  usual  Kalman  filter  using  guaranteed  symmetry  of  the  covariance 
matrix.  Hendrickson  reports  steady  state  nf.vigation  accuracies  wh;.  h  vary 
between  1000  ft  in  the  vicinity  of  the  landmark  to  3000  ft  during  landmark 
voids . 

Because  of  the  comprehensiveness  of  the  above  mentioned  studies,  no  attempt 
was  made  in  the  current  effort  to  analyze  unknown  landmark  tracking  unless  it 
is  part  of  a  hybrid  system. 

However,  as  a  check  on  the  SNAP  software,  several  ONAP  simulations  identical 
to  those  reported  in  the  references  were  performed.  Comparison  of  the  results 
showed  good  agreement. 
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5.  HORIZON  SENSORS 


5.  1  DISCUSSION 

Autonomous  navigation  with  horizon  sensors  appears  to  be  attractive  for  the 
following  reasons: 

a.  Complete  Autonomy 

Other  than  a  knowledge  of  inertial  attitude  and  time,  no  other 
information  is  required  for  horizon  sensors  to  provide  complete 
navigation  information  in  three  dimensions. 

b.  Low  Cost 

Compared  to  other  space  navigation  sensors  and/or  systems, 
the  cost  of  horizon  sensors  is  quite  attractive  since  the  total 
system  only  requires  horizon  sensors,  a  clock,  an  attitude  ref¬ 
erence  system,  and  a  computer  to  process  the  measurements. 

c.  Deterministic  Operation 

If  properly  used,  horizon  sensors  can  yield  a  deterministic  posi¬ 
tion  fix  in  three  dimensions.  Such  information  could  be  valuable 
for  bootstrap  initialization  of  more  sophisticated  data  processing 
algorithms.  This  property  can  also  be  useful  after  dormant 
navigation  periods  or  for  recovery  from  possible  spacecraft 
maneuvers. 

d.  Wide-Range  Linear  Operation 

Recovery  from  large  initial  error  is  possible  with  straightforward 
application  of  usual  filter  algorithms  because  large  position  dis¬ 
placements  result  in  only  small  measurement  angle  changes. 

Unfortunately,  these  advantages  are  offset  by  the  uncertainty  in  the 
horizon  definition  at  any  particular  time.  This  uncertainty  is  a  dominant 
error  source  during  low-altitude  operation.  Currently,  the  lack  of  an  ade¬ 
quate  model  for  this  error  source  is  a  significant  limiting  factor  in  horizon 
sensor  navigation  performance  during  low  orbit  operation. 


The  use  of  horizon  sensors  for  navigation  has  received  considerable  attention 
by  others  (Refs.  15,  20,  24,  28,  35,  and  46).  The  most  commonly  used  infor¬ 
mation  provided  by  horizon  sensors  is  the  measurement  of  the  direction  to  the 
center  of  the  earth.  Since  this  measurement  by  itself  only  provides  position 
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information  in  the  local  horizontal  plane,  altitude  can  be  deduced  by  processing 
a  time  sequence  of  measurements  with  the  known  dynamical  model.  As  a  sec¬ 
ond  alternative,  a  radar  altimeter  can  be  used  to  directly  measure  the  altitude. 
This  may  not  be  desirable  for  autonomous  navigation  if  it  is  preferable  that 
the  navigator  be  passive.  A  third  possible  choice  for  determining  altitude 
involves  the  use  of  information  inherently  present  in  horizon  sensor  measure¬ 
ments.  Specifically,  the  angle  subtended  by  the  widest  apparent  diameter  of 
the  earth's  disk  is  readily  computed  from  horizon  sensor  data  and  provides  a 
direct  approach  for  determining  altitude.  This  is  referred  to  as  a  planet 
diameter  measurement  in  Ref.  3.  If  it  were  necessary,  either  the  radar 
altimeter  or  the  angle  subtended  by  the  earth's  disk  could  provide  determinis¬ 
tic  position  information  in  three  dimensions  when  combined  with  the  measure¬ 
ment  of  the  direction  to  the  center  of  the  earth.  Although  a  precise  knowledge 
of  the  sensor  orientation  in  inertial  space  is  needed,  no  a  priori  position 
information  is  strictly  required. 

The  geometry  of  horizon  sensor  measurements  is  shown  in  Figure  5-1.  The 
right-handed  cartesian  coordinate  system  in  Figure  5-1  is  oriented  such  that 
the  z2-axis  nominally  points  opposite  to  the  direction  of  local  vertical  and  Zj 
and  z3  are  in  the  local  horizontal  plane.  The  inability  of  to  point  exactly 
opposite  to  the  local  vertical  direction  is  limited  by  the  accuracy  of  the  atti¬ 
tude  control  system  of  the  vehicle  and  has  no  particular  bearing  on  the  navi¬ 
gation  problem,  providing  the  errors  are  not  excessively  large.  Although  the 
z-coordinate  system  is  only  coarsely  aligned  with  the  local  vertical,  its  atti¬ 
tude  in  inertial  space  is  assumed  to  be  known  by  some  precise  attitude  deter¬ 
mination  system  such  as  SPARS  (Refs.  10  and  11).  Hence,  the  attitude  errors 
are  assumed  to  be  negligible  compared  to  the  errors  in  the  horizon  sensor. 

The  measurement  of  the  direction  to  the  center  of  the  earth  is  defined 
by  two  angles  0^  and  0j  which  are  shown  in  Figure  5-2a.  These  angles  can 
by  expressed  in  terms  of  the  ECI  (see  Appendix  C)  vehicle  position  coordinates, 


5-2 


Figure  5-2a.  Horizon  Sensor  Measurement  Geometry 


(5-2) 


0'j  =  arc  cos 
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and  the  rotation  matrix  B  is  known  from  the  assumed  attitude  reference 
system.  Figure  5- 2b  has  been  drawn  to  emphasize  the  angle  definitions.  In 
reality,  0j  and  0J  should  be  very  close  to  90  deg  if  the  attitude  control  system 
is  performing  properly.  The  angle  subtended  by  the  widest  diameter  of  the 
earth's  apparent  disk  is  defined  to  be  @2  and  is  shown  in  Figure  5-2b;  can 
be  expressed  as 


©2=2  arc  sin  (Re/(Re  +  h)]  (5-4) 

where  R  is  the  radius  of  the  earth  and  h  is  the  altitude  of  the  navigator. 

These  measurements  are  corrupted  by  both  sensor  noise  and  horizon  definition 
errors.  The  former  is  due  to  the  inability  of  the  sensor  to  perfectly  discern 
the  14  -  16  |jl  CO^  absorption  layer  (the  phenomena  detected  by  high  perfor¬ 
mance  sensors  such  as  Quantic  Mod  IV)  (Ref.  46).  The  magnitude  of  this 
error  is  on  the  order  of  0.01  deg  (l<r)  and  is  assumed  to  be  white  noise  in 
this  study.  The  dominant  error  source,  however,  is  not  in  the  sensors 
themselves,  but  in  the  uncertainty  of  the  altitude  of  the  CO2  absorption  layer. 

It  can  be  shown  that  these  horizon  errors  are  time-  and  space-correlated 
(Refs.  16  and  34).  Thus  they  can  be  modeled  as  Markov  processes  (Ref.  28). 
The  derivation  of  the  appropriate  relationships  describing  the  effect  of  these 
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errors  on  the  observations  is  somewhat  involved  and  only  the  key  results  are 
given  here.  It  is  shown  in  Appendix  D  that  the  horizon  errors  in  0j,  0  j ,  and 
0£,  represented  by  the  state  variables  x^,  Xg,  Xg,  respectively,  propagate  in 
time  according  to  the  differential  equations 

x^  =  -x.y/625  +  u7 

xg  =  -Xg/625  +  Ug  (5-5) 

Xg  =  "Xg/625  +  Ug 

where 

£[x*(0)]  =  0.  014852  deg2 

£[x2(0)]  =  0.014852  deg2  (5-6) 

£[x2(0)]  =  0.051452  deg2 

and  Uy,  Ug,  and  Ug  are  white  noise  defined  by 

£[u7(t)u7(r)]  =  7.  05  X  i 0”?6(t  -  t)  deg2/sec2 

£[ug(t)ug(*r)]  =  7.  05  X  1 0_76(t  -  t)  deg2/ sec2  (5-7) 

^[ug(t)ug(T)]  =  8.  47  X  1 0~^6(t  -  t)  deg2/ sec2 

5.  2  BASELINE  CONFIGURATION  AMD  .ASSOCIATED  RESULTS 

The  system  discussed  in  Sect.  5.  1  was  analyzed  in  a  series  of  simulations. 
The  baseline  configuration  is  described  in  Table  5-1. 
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Table  5-1.  Baseline  Horizon  Sensor  Configuration 


Key  Parameters 


Measurement  Model 

1.  Angle  subtended  by  disk  of  earth 

2.  Two  angles  to  apparent  center  of  earth 

Gravity  Model 

Real  World:  APL  "8-8"  ! 

Filter  World:  p  and  ^ 

Atmosphere  Model 

Real  World:  U.S.  Standard  1962 

Filter  World:  Exponential 

Data  Rate 

1  Sample/30  sec 

Measurement  Noise 
and  Biases 

0.  01  deg  { 1  tr )  random  noise 

Markov  biases  (see  Sect.  5.  1) 

Forcing  Noise 

None  except  for  noise  on  Markov  processes 
(see  Eqs.  (5-5),  (5-6),  (5-7)) 

Initial  Conditions 


States 

Real  World 

Filter  World 
(all  states 
estimated) 

Standard 

Deviation 

Xj,X2,  Xj  (inertial 
position) 

Reference 

Orbit 

6700  ft  offset 
(in  each  direction) 

10000  ft 

x.,,xe,Xr  (inertial 

4  5  6  velocity) 

Reference 

Orbit 

6.7  ft/ sec  offset 
(in  each  direction) 

1 0  ft/sec 

x_,x„  (Markov 
biases) 

0  deg 

0.  01485  deg 

0.  01485  deg 

Xg  (Markov  bias) 

0  deg 

0. 05145  deg 

0. 05145  deg 
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5.  2.  1  The  results  of  simulating  the  baseline  configuration  are  shown 
graphically  in  Figure  5-3  and  are  indicated  below.  After  about  one-half 
orbit  the  steady  state  condition  is  reached.  The  position  standard  deviations 
converge  to 

1400  ft  Radial 
4500  ft  Intrack 
2800  ft  Crosstrack 
5500  ft  RSS 

5.  2.  2  Model  errors  do  not  appear  to  have  an  effect  until  the  second  orbit. 

5.2.  3  The  Markov  biases  in  the  sensor  measurements  are  weakly  observable. 
The  standard  deviations  of  the  bias  in  the  angles  defining  the  direction  to  the 
center  of  the  earth  were  reduced  to  2/3  of  their  a  priori  values  after  half  an 
orbit,  while  the  standard  deviation  for  the  bias  in  the  angle  subtended  by  the 
earth's  disk  was  reduced  by  1/2  in  that  same  time  interval. 

5.2.4  The  basic  performance  of  the  horizon  sensor  is  limited  by  the  error 
model  used  for  the  altitude  of  the  CO^  absorption  layer. 

5.  3  INITIAL  CONVERGENCE  ANALYSIS 

5.3.1  The  baseline  case  was  compared  to  cases  with  initial  position  errors 
and  corresponding  standard  deviations  of  5,  10,  and  20  miles.  Table  5-2  gives 
the  position  errors  obtained  after  the  angle  subtended  by  the  earth  disk  and  two 
angles  defining  the  direction  to  the  center  of  the  earth  have  been  processed  by 
the  filter.  The  agreement  between  the  state  residuals  and  the  corresponding 
standard  deviations  indicates  that  a  linear  algorithm  is  sufficient  for  initiali¬ 
zation  for  position  offset  as  large  as  20  miles.  With  errors  this  large,  the 
effect  of  the  a  priori  covariance  is  negligible.  Thus  a  deterministic  solution 
should  give  comparable  errors. 
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Figure  5-3.  RSS  Position  Errors  for  Baseline  Horizon  Sensor  System 


Table  5-2.  Errors  Obtained  after  One  Set  of  Horizon  Sensor  Measurements 


RSS 

Standard 

Deviation 

(ft) 

8607 

10212 

10311 

10394 

State 

Residual 

(ft) 

6882. 

8496 

7762 

6752 

Crosstrack 

Standard 

Deviation 

(ft) 

5600 

6800 

6900 

6900 

State 

Residual 

(ft) 

•  •  •  • 

+  in  f* 

•  ^  O' 

^  W  O' 

sO  o  m 

i  i  i  i 

Intrack 

Standard 

Deviation 

(ft) 

5600 

6800 

6900 

6900 

State 

Residual 

(ft) 

-974. 

-5104. 

-4717. 

-3960. 

*■4 

■H 

•o 

Standard 

Deviation 

(ft) 

3100 

3300 

3400 

3500 

It 

05 

State 

Residual 

(ft) 

-3004. 

-2500. 

-1625. 

1699. 

Initial  Position 
Errors  (n  mi) 

1 

5 

10 

20 

5 


5.  3.  2  Convergence  to  the  values  given  in  Sect.  5.  2. 1  occurs  in  about 
1/2  orbit  and  is  apparently  independent  of  the  magnitude  of  the  initial  errors 
for  the  particular  orbit  studied. 

5.3.3  Removal  of  the  information  provided  by  the  measurement  of  the  angle 
subtended  by  the  earth's  disk  resulted  in  a  delay  of  convergence  until  the  end 
of  the  first  orbit.  Furthermore,  such  a  system  does  not  provide  the  capa¬ 
bility  for  deterministic  operation. 

5.3.4  Although  the  accuracy  attainable  with  horizon  sensors  is  limited,  the 
advantages  gained  by  their  convergence  properties  suggest  their  use  as  auxil¬ 
iary  sensors  in  a  general  autonomous  navigation  system. 
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6.  UNKNOWN  LANDMARK  TRACKER  WITH  HORIZON  SENSORS 


Upon  examining  the  results  on  both  unknown  landmark  tracking  and  horizon 
sensors,  it  is  natural  to  consider  a  hybrid  configuration  containing  both 
sensors.  The  properties  of  these  sensors  certainly  appear  to  complement 
each  other  quite  well.  Accurate  navigation  information  could  be  supplied  by 
the  unknown  landmark  tracker  while  initialization  and  error  bounding  during 
landmark  voids  could  easily  be  accommodated  by  the  horizon  sensor.  In 
addition,  both  sensors  are  highly  autonomous. 

6.  1  BASELINE  CONFIGURATION  AND  ASSOCIATED 

RESULTS 

The  baseline  configuration  used  is  described  in  detail  in  Table  6-1. 

6.1.1  The  results  obtained  by  simulating  the  baseline  case  are  indicated  graph¬ 
ically  in  Figure  6-1  and  summarized  in  Table  6-2.  Although  the  performance 
is  about  three  times  better  than  with  the  horizon  sensor  alone,  the  system 
does  not  appear  capable  of  producing  accuracies  much  better  than  2000  ft. 

6.  1.2  Removal  of  the  horizon  sensors  from  the  baseline  configuration  yields 
the  results  shown  in  Figure  6-2  and  Table  6-3.  Comparison  of  Tables  6-2 
and  6-3  shows  how  valuable  horizon  sensors  can  be  for  initialization  and  for 
bounding  error  growth  during  landmark  voids. 

6.  2  EFFECT  OF  INCREASED  LANDMARK 

VISIBILITY  TIME 

As  a  possible  improvement  to  the  baseline  performance,  the  line  of  sight- 
nadir  angle  was  increased  from  45  to  60  deg  to  allow  longer  visibility  times 
for  each  landmark.  Atmospheric  refraction  was  ignored  to  simplify  the 
analysis.  A  comparison  of  the  results  displayed  in  Table  6-4  with  those  of 
Table  6-2  shows  the  effect  is  relatively  marginal.  Indeed,  since  the  hard¬ 
ware  problem  associated  with  landmark  definition  near  the  horizon  is  quite 
formidable,  there  is  little  reason  to  consider  the  approach  any  further. 
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Table  6-1 


Baseline  Unknown  Landmark  Tracker-Horizon 
Sensor  Configuration 


Kay  Parameters 


Meaaurement  Set 

1.  Two  angles  to  an  unknown  landmark. 

2.  Two  angles  to  the  apparent  center  of  the  earth 
(biaecto  angles). 

3.  Angle  subtended  by  earth's  disk 

Gravity  Model 

Real  World:  APL  8-8 

Filter  World:  pend  Jj 

Atmosphere  Model 

Real  World:  U.S.  Standard  1962 

Filter  World:  exponential 

Data  Rate 

Unknown  landmark  tracker:  Landmarks  were  "generated" 
in  the  orbit  plane  45  deg  from  nadir  at  120,  450,  1140, 

1620.  3360,  and  3900  sec  from  epoch.  Observations 
were  taken  every  10  sec  until  the  line  of  sight  to  the 
landmark  was  more  than  45  deg  from  nadir. 

Horizon  sensors:  1  sample  30/sec  (continuously) 

Meaaurement  Noise  and  Biases 

Unknown  landmark  tracker:  10  arc  sec  random  (la), 

30  arc  sec  bias 

Horizon  sensor:  0.01  degree  random  (la) 

Markov  biases  (see  Sect.  5.  1) 

Forcing  Noise 

None  except  for  Markov  biases. 

Initial  Condition* 


States 

Real  World 

Filter  World 
(all  states  estimated) 

Standard 

Deviation 

Xj,  Xj,  x.j  (Inertial  position) 

Reference  Orbit 

6700  ft  offset 
in  each  direction 

10,000  ft 

x^,  Xg,  x^  (Inertial  velocity) 

Reference  Orbit 

6.  7  ft/sec  offset 
in  each  direction 

10  ft /sec 

x?,  Xg  (Markov  biases) 

0  deg 

0, 01485  deg 

0.  01485  deg 

Xg  (Markov  bias) 

0  deg 

0,  05145  deg 

0.  05145  deg 

x10>  Xjj  (Landmark  latitude,  longitude) 

Computed 

Computed 

Computed 

Xj2  (Landmark  altitude) 

3000  ft 

Random  offset 

1500  ft  (la) 

1500  ft  (la) 

x13'  x14  angle  biases: 

0  deg 

0.  0083  deg 

0,  0083  deg 

fi 


(09J  JO  spuosnom )  yoad3  NOlilSOd  SSH 


6-3 


Figure  6-i.  RSS  Position  Errors  for  Baseline  Unknown  Landmark 
Tracker  -  Horizon  Sensor  System 


Table  6-2.  RSS  Position  Error  for  Baseline  Unknown 
Landmark  Tracker-Horizon  Sensor 
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‘igure  6-2.  RSS  Position  Errors  for  Unknown  Landmark  Tracker 


Table  6-3.  RSS  Position  Error  for  Baseline  Unknown  Landmark 
Tracker  without  Horizon  Sensor 


Time  of  Acquisition 
and  Termination  of 
Landmark  (sec) 


Landmark 

Number 


RSS 

Error 

(ft) 

RSS 

Standard 

Deviation 

(ft) 

11,640 

17,321 

7,851 

2,392 


7, 446 
3,  888 


4,201 
3,  503 


6,  291 
2,732 


24,814 

2,911 


3,53 

2,69 


2,647 

1,817 


5400 


3,  178 
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Table  6-4.  Effect  of  60  Deg  Maximum  Line  of  Sight-Nadir 
Angle  on  Unknown  Landmark  Tracker-Horizon 
Senior  RSS  Position  Error 


RSS 

Time  of  Acquisition 

RSS 

Standard 

and  Termination  of 

Landmark 

Error 

Deviation 

Landmark  (sec) 

Number 

(ft) 

(ft) 

0 

- 

11.  604 

17,321 

120 

1 

n 

7,  570 

210* 

HI 

6,  766 

450 

5,  884 

540* 

l 

3,  853 

1140 

2,950 

2,823 

1200* 

J 

1,322 

2,  117 

1620 

A 

1,567 

2,225 

1680* 

706 

1,903 

3360 

e 

2,429 

3,  650 

3450* 

D 

1,661 

1,719 

3900 

A 

1,591 

1,538 

4020* 

2,092 

1,285 

5400 

- 

1,459 

2,  082 

These  entries  are  not  landmark  termination  times. 
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EFFECT  OF  LANDMARK  DENSITY 


6.  3 

6.  3.  1  In  order  to  determine  the  effect  of  increasing  the  landmark  denaity, 
several  two-orbit  computer  runs  were  made.  First,  the  baseline  simulation 
was  extended  to  the  end  of  the  second  orbit;  then,  a  similar  configuration 
with  twelve  more  or  less  uniformly  distributed  landmarks  per  orbit  was 
simulated;  and,  finally,  the  hypothetical  situation  of  unlimited  landmark 
availability  was  considered.  In  all  cases,  the  actual  RSS  error  became 
meaningless  in  the  second  orbit  because  of  uncompensated  model  errors. 

Thus  only  the  RSS  standard  deviations  were  reported.  The  results,  displayed 
in  Table  6-5,  represent  a  lower  bound  since  the  effect  of  model  error  is  not 
included. 

6.  3.2  Basically,  the  results  show  that  a  major  increase  in  landmark  density 
is  required  to  bring  about  a  significant  decrease  in  error.  It  is  doubtful  that 
a  sufficient  increase  in  landmarks  could  be  obtained  in  practical  situations  to 
warrant  the  effort. 

6.4  EFFECT  OF  LANDMARK  DISTRIBUTION 

The  baseline  landmark  acquisition  schedule  chosen  did  not  allow  for  possible 
voids  due  to  operation  over  cloud-covered  or  night  areas.  In  Refs.  5,  12, 
and  16,  the  day-night  problem  was  considered  using  five  landmarks  per 
orbit,  all  concentrated  on  the  "day"  side.  The  results  reported  showed 
oscillations  in  accuracy  from  1000  to  3000  ft.  As  shown  in  Table  6-5,  use 
of  a  similar  number  of  landmarks  distributed  more  evenly  prevents  such 
large  variations  with  very  little  degradation  in  maximum  attainable  accuracy. 
Thus,  landmark  distribution  is  critical  for  bounding  errors,  but  is  not  very 
effective  for  improving  the  navigation  accuracy. 

6.  5  USE  OF  LANDMARKS  WITH  KNOWN  ALTITUDES 

6.  5.  1  If  landmarks  were  restricted  to  the  coastline,  the  uncertainty  in 
landmark  altitude  would  become  extremely  small.  The  effect  is  essentially 
equivalent  to  the  removal  of  landmark  altitude  from  the  state  vector.  The 
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Table  6-5.  Effect  of  Landmark  Density  on  Unknown  Landmark 
Tracker-Horizon  Sensor  RSS  Position  Error 


results  obtained  from  simulating  such  a  situation  are  sumarized  In  Table  6-6. 
Comparison  with  Table  6-2  shows  that  a  SO  percent  improvement  can  be 
obtained  by  the  suggested  technique. 

6.  5.  2  As  might  be  expected,  the  use  of  landmarks  with  known  altitudes 
reduces  the  radial  error  drastically.  An  order  of  magnitude  reduction 
in  this  error  was  observed. 

6.6  PERFORMANCE  WITH  A  HIGH  ACCURACY 

SENSOR 

6.  6.  1  In  order  to  assess  the  "ultimate"  performance  of  an  unknown  land¬ 
mark  tracker,  a  high  accuracy  sensor  with  2  arc  sec  of  random  noise  (l<r) 
and  4  arc  sec  of  bias  was  considered.  The  results  of  using  this  configura¬ 
tion  are  shown  in  Column  A  of  Table  6-7.  While  the  errors  are  considerably 
less  than  those  observed  for  the  baseline  case,  they  were  not  significantly 
different  from  those  using  the  baseline  sensor  with  coastline  landmarks. 

6.  6.  2  The  effect  of  assuming  coastline  landmark  with  the  high  accuracy 
sensors  is  shown  in  Column  B  of  Table  6-7.  A  considerable  improvement 
in  performance  of  the  error  analysis  is  clearly  obtained.  Note,  however, 
that  model  errors  now  become  quite  significant  by  the  end  of  the  first  orbit. 

6.  6.  3  The  effect  of  removal  of  both  atmosphere  and  geopotential  model 
errors  is  shown  in  Column  C  of  Table  6-7.  These  results  indicate  that 
improved  geopotential  models  and  the  possible  use  of  accelerometers  for 
drag  measurement  could  provide  a  significant  increase  in  accuracy. 

6.  6.4  It  should  be  noted  that  although  high  accuracy  navigation  with  unknown 
landmark  trackers  is  indicated  by  these  results,  the  underlying  assumptions 
are  not  realistic.  The  feasibility  of  building  a  sensor  with  the  required 
accuracy  and  the  ability  to  distinguish  coastlines  from  other  landmarks  has, 
by  no  means,  been  demonstrated.  Furthermore,  the  results  obtained  were 
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Table  6-6.  RSS  Position  Errors  for  Baseline  Unknown  Landmark 
Tracker  and  Horizon  Sensor  with  Coastline  Landmarks 


RSS 

Time  of  Acquisition 

RSS 

Standard 

and  Termination  of 

Landmark 

Error 

Deviation 

Landmark  (sec) 

Number 

(ft) 

(ft) 

0 

. 

11,640 

1732 

120 

1 

13,230 

7570 

210 

14,  646 

6638 

450 

2 

13,  677 

5834 

540 

5,475 

2856 

1140 

i 

3,  729 

2738 

1200 

o 

1.  530 

1320 

1620 

A 

119 

1095 

1680 

* 

328 

974 

3360 

c 

2,  574 

1349 

3450 

D 

1,937 

1189 

3900 

L 

1.  192 

1047 

4020 

0 

763 

914  | 

5400 

- 

2,  061 

964 
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4020  592  613  637  185  220 


baaed  on  a  relatively  uniform  landmark  distribution,  whereas  a  real  system 
must  allow  for  sizable  landmark  voids  during  operation  at  night  or  under 
conditions  with  heavy  cloud  cover. 
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7.  RANGE  RATE  TO  KNOWN  GROUND  STATIONS 


Although  the  configuration  considered  in  the  previous  section  was  highly 
autonomous,  certain  rather  optimistic  assumptions  were  required  to  obtain 
high  accuracy  navigation  performance.  In  an  attempt  to  meet  such  a  goal 
more  realistically,  the  autonomy  requirement  was  relaxed  and  navigation  by 
means  of  range  rate  to  known  ground  stations  was  considered. 

BASELINE  CONFIGURATION  AND  ASSOCIATED 


RESULTS 


The  baseline  configuration  described  in  Table  7-1  used  21  ground  stations 
which  were  chosen  at  locations  near  existing  doppler  installations.  The 
assumed  coordinates  are  given  in  Appendix  F. 


7.1.1  The  ground  trace  of  the  first  orbit  for  the  reference  ephemeris  of 
Appendix  C  is  shown  in  Fig.  7-1.  Those  ground  stations  visible  to  the  satellite 
are  also  indicated  on  that  diagram. 

7.1.2  The  time  history  of  the  RSS  position  errors  for  the  baseline  case  is  shown 
graphically  in  Figure  7-2.  The  results  before  and  after  acquisition  of  each 
ground  station  are  summarized  in  Table  7-2.  Although  the  errors  decrease 
after  each  ground  station  acquisition,  they  grow  quite  rapidly  when  no 
measurements  are  taken.  By  the  second  orbit,  however,  this  pattern  of 

wild  error  fluctuation  appears  to  stabilize. 

7.1.3  Unmodeled  dynamic  errors  do  not  become  significant  until  well  into  the 
second  orbit. 

7.2  EFFECT  OF  INCREASING  NUMBER  OF 

GROUND  STATIONS 

7.2.1  It  was  hypothesized  that  the  relative  instability  noted  in  the  first  orbit  of 
the  baseline  case  can  be  attributed  to  an  insufficient  number  of  ground  station 
acquisitions.  In  order  to  test  this  hypothesis,  additional  stations  were 
assumed.  In  choosing  the  locations,  an  attempt  was  made  to  provide  as 
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Table  7-1.  Baseline  Configuration  for  Range  Rate  to  Known  Ground  Stations 


Key  Parameters 


Measurement  Model 

Instantaneous  range  rate 

Number  of  Stations 

21 

Gravity  Mo^el 

Real  World:  8-8 

Filter  World:  jjl  and 

Atmosphere  Model 

Real  World:  U.S.  Standard  1962 
Filter  World:  exponential 

Data  Rate 

1  sample/10  sec 

Measurement  Noise 

1 . 67  ft/sec  (Id) 

Forcing  Noise 

None 

Minimum  Elevation  Angle 
for  ground  station  visibility 

5  deg 

Initial  Conditions 


State 

Real  World 

Filter  World 
(all  states  estimated) 

Standard  Deviation 

Xj.X3.X3 

(inertial  position) 

Reference  Orbit 

6700  ft  offset 

(in  each  direction) 

10,  000  ft/sec 

x4.x5.x6 

(inertial  velocity) 

Reference  Orbit 

6.  7  ft/sec  offset 

(in  each  direction) 

10  ft/ sec 

7-2 


i 


*70  110  170  1*0  1*0  140  110  IIP  110  100  *0  10  70  *0  SO  40  *0  10  10  0 

LONGITUDE 


Figure  7-1.  First  Orbit  Ground  Trace  Showing  Visible  Stations 
(21  ground  station  system' 
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Figure  7-2.  RSS  Position  Errors  for  Baseline  Range-Rate  System 


Table  7 -2.  RSS  Position  Errors  for  Baseline 
Range  Rate  Configuration 


Time  of  Acquisition 
and  Termination  of 
Ground  Station  (sec) 


Station 

Number 


RSS 

Error  (ft) 


11605 


RSS 

Standard 

Deviation 

(ft) 


17321 


28938 


near  worldwide  coverage  as  possible,  without,  of  course,  utilizing  any 
politically  excluded  territory.  Island  installations  were  used  extensively 
to  avoid  large  gaps  in  ocean  coverage.  The  final  list  contains  104  stations, 
including  the  21  used  for  the  baseline  case  (see  Appendix  F).  The  coordinates 
of  stations  22  through  104  were  rounded  to  the  nearest  5  deg  to  facilitate 
program  input  and  graphical  work.  The  first  orbit  ground  trace  with  the 
associated  visible  stations  using  the  complete  list  is  shown  in  Figure  7-3. 

7.  2.  2  The  results  of  a  range  rate  navigation  simulation  with  the  104  station 
deployment  are  displayed  graphically  in  Figure  7-4.  As  anticipated,  error 
osillations  are  drastically  reduced  and  a  low  steady  state  error  is  reached 
quite  rapidly  in  about  one- third  of  an  orbit.  However,  because  of  the  improved 
accuracy,  model  errors  become  noticeable  during  the  first  orbit. 

7.  3  INITIAL  CONVEF  GENCE  ANALYSIS 

Cases  with  initial  position  offsets  of  5  and  10  n  mi  with  104  ground  stations 
were  studied  in  order  to  determine  the  convergence  properties  of  the  assumed 
measurement  configuration.  The  results,  which  are  summarized  in  Table  7-3, 
show  that  filter  convergence  with  initial  errors  of  5  n  mi  or  more  cannot  be 
obtained.  This  suggests  possible  augmentation  of  the  range  rate  configuration 
with  horizon  sensors. 

7.  4  RANGE  RATE  WITH  ANGLE  SUBTENDED  BY 

EARTH'S  DISK 

7.4.  1  Since  range  rate  observations  do  not  require  precise  vehicle  attitude 
information,  it  would  be  desirable  that  any  sensor  introduced  for  improving 
convergence  properties  should  also  be  independent  of  such  information.  While 
the  complete  horizon  sensor  system  described  in  Section  5  implicitly  requires 
a  precise  attitude  reference  system,  the  angle  subtended  by  the  earth's  disk 
does  not,  and  could  conceivably  be  used  to  augment  the  range  rate  measurement. 
This  contention  was  tested  in  a  series  of  simulations  in  which  the  initial  posi¬ 
tion  offset  was  varied.  The  results  are  summarized  in  Table  7-4. 

7.4.2  Comparison  of  Tables  7-4  and  7-2  shows  that,  while  the  angle  subtended 
by  the  earth's  disk  prevents  catastrophic  failures,  it  by  no  means  provides  an 
adequate  solution  to  the  convergence  problem. 
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Figure  7-4.  RSS  Position  Errors  Range-Rate  System  with  104  Ground  Stations 


Table  7-3.  RSS  Position  Errors  for  Varicus  Initial  Offsets  Using  Baseline  Range  Data 
Configuration  With  104  Stations 


3 
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Table  7-4.  RSS  Position  Errors  for  Various  Initial  Offsets  Using  Range  Rate 
Earth  Disk  System  With  104  Stations 
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RANGE  RATE  WITH  COMPLETE  HORIZON 


SENSOR  SYSTEM 

From  the  above  results,  it  appears  that  the  full  horizon  sensor  system  of 
Section  5  should  be  investigated  as  a  means  of  preventing  convergence 
problems.  The  results  of  combining  the  baseline  horizon  sensor  and  range 
rate  systems  are  summarized  in  Table  7-5.  By  comparing  the  results  with 
those  of  Table  7-2,  the  utility  of  the  horizon  sensor  as  an  effective  means  of 
bounding  errors  is  quite  evident.  Furthermore,  no  significant  degradation 
in  behavior  was  noted  when  the  system  was  faced  with  large  initial  errors. 
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Table  7-5.  RSS  Position  Errors  With  Range  Rate  and  Horizon 
Sensor  Measurements 


Time  of  Acquisition 

RSS 

RSS  Standard 

and  Termination  of 

Error 

Deviation 

Ground  Station  (sec) 

Station  Number 

(ft) 

(ft) 

0 

-- 

11605 

17321 

930 

19 

11999 

6963 

990 

2000 

3240 

1170 

4 

2010 

3326 

1320 

782 

2003 

1560 

3 

2843 

2394 

1800 

2320 

1327 

3180 

10 

12761 

3398 

3390 

263 

457 

5400 

-- 

1536 

2359 

8.  ONE-WAY  DOPPLER  MEASUREMENTS  TO  KNOWN 

GROUND  STATIONS 


8.  1  DISCUSSION 

As  shown  in  Sect.  7,  range  rate  measurements  to  known  landmarks  appear 
to  provide  extremely  accurate  satellite  navigation.  Actually,  range  rate  is 
a  mathematical  abstraction  of  a  real  process.  If  the  results  obtained  thus  far 
are  to  be  meaningful,  that  process  must  be  examined  carefully. 

Specifically,  range  rate  data  are  actually  data  derived  from  doppler  shifts  of 
high  frequency  electromagnetic  radiation.  Practical  systems,  yielding  what 
are  commonly  termed  "one-way  doppler  measurements",  consist  of  several 
ground  stations  located  around  the  globe  which  continuously  transmit  a  nominal 
carrier  frequency  (Refs.  45  and  52).  When  one  of  the  ground  stations  becomes 
visible  to  a  user  satellite,  the  frequency  of  the  signal  received  by  the  satellite 
is  compared  with  the  frequency  of  a  stable  onboard  oscillator.  The  compari¬ 
son  ideally  yields  a  difference  frequency  Af  that  is  directly  proportional  to 
the  time-rate  of  change  of  the  distance  between  the  satellite  and  the  ground 
station,  (p),  that  is, 


p  =  kAf 


(8-1) 


However,  since  frequency  is  defined  as  the  number  of  events  per  unit  time, 
measurement  of  the  desired  frequency  difference  cannot  be  made  instanta¬ 
neously.  Instead,  each  measurement  must  be  taken  over  some  finite  interval 
t.  S  t  S  tj.  In  practice,  the  integral  of  the  frequency  difference  over  that 
interval  is  measured.  Mathematically,  the  observation  can  be  represented  as 


(8-2) 
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Although  Eq.  (8-2)  can  be  divided  by  -  t^  to  yield  the  average  range  rate 
over  that  interval,  it  is  more  convenient  to  consider  the  range  difference  of 
Eq.  (8-2)  directly.  In  any  case,  it  should  be  clear  that  instantaneous  range 
rate  is  only  an  approximation  to  the  actual  data  type. 

In  actual  practice  the  satellite  and  ground  frequency  oscillators  are  not  per¬ 
fect.  Over  short  intervals  the  error  can  be  modeled  with  a  range  bias,  a 
frequency  offset  and  a  frequency  offset  drift  with  white  noise  assumed  on 
both  the  range  bias  and  the  frequency  offset  (Ref.  41).  Thus,  if  the  state 
vector  is  defined  as 


x2 

x3 


EC1  satellite  position 


x4 

x5 

x6 


ECI  satellite  velocity 


instantaneous  range  bias,  (ft) 
instantaneous  frequency  offset,  (ft/ sec) 


ECI  ground  station  position 


x  1 2 
X13 
X 14 


ECI  ground  station  velocity 
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x j g  range  to  satellite  at  acquisition  time  (see  Appendix  E) 

/  2 

rate  of  change  of  offset  in  frequency  (ft/ sec  ) 


where 


*7  =  *8  +  u 
*8  *  *16 
*16  =  0 


and  u  is  white  noise, 


the  measurement  can  be  expressed  as 


y  =  p(t{)  +x?(tf)  -  p(tt )  -  x?(t.)  +  v(t£)  -  v(t.) 


(8-3) 


(8-4) 


where  v  is  white  noise.  The  error  model  of  Eq.  (8-3)  may  be  represented 
as  shown  in  Fig.  8-1. 


Figure  8-1.  Block  Diagram  for  One-Way  Doppler  Error  Model 

The  measurements  of  Eq.  (8-4)  represent  an  additional  problem  when  Kalman 
filtering  processing  is  required.  As  formulated,  the  Kalman  filter  algorithm 
(Refs.  18,  19,  20,  and  29)  requires  that  the  observations  be  expressed  as  a 


function  of  the  state  vector  at  a  specific  time.  The  observation  described  in 
Eq.  (8-4),  however,  involves  the  state  at  two  different  times,  i.e.,  t.  and  t^. 

Thus  an  appropriate  technique  must  be  used  to  modify  Eq.  (8-4)  so  that  com¬ 
patibility  with  the  Kalman  filter  can  be  obtained. 

Three  approaches  were  considered  to  resolve  this  difficulty.  The  first  makes 
use  of  the  range  at  acquisition  as  a  state  variable  (Acquisition  Range  Reference 
Method),  the  second  uses  state  extrapolation  and  is  based  on  the  difference 
between  ranges  at  two  successive  observation  times  (Previous  Range  Refer¬ 
ence  Method),  while  the  last  method  is  based  on  range  rate  as  an  approxi¬ 
mation  to  the  actual  measurement  (Range  Rate  at  Midpoint).  A  detailed 
development  of  these  methods  is  provided  in  Appendix  E. 

8.  2  BASELINE  CONFIGURATION  AND  ASSOCIATED  RESULTS 

The  baseline  configuration  used  to  examine  one-way  doppler  systems  is 
described  in  detail  in  Table  8-1. 

8.  2.  1  The  time  history  of  the  RSS  position  errors  for  the  baseline  case  is 
shown  graphically  in  Figure  8^2.  A  similar  plot  is  given  in  Figure  8-3  for 
a  case  with  the  complete  104  ground  stations  of  Appendix  F. 

8.  2.  2  The  results  obtained  are  very  similar  to  the  range-rate  results  of  Sect.  7. 

8.  2.  3  After  1-1/3  orbits  with  the  baseline  case,  position  errors  are  reduced  to 
about  700  ft.  After  that,  model  errors  become  significant. 

8.  2.  4  After  1/3  orbit,  with  the  104  station  case,  the  position  errors  are 
reduced  to  about  300  ft.  Thereafter,  model  errors  become  noticeable. 

8.  2.  5  As  shown  in  Appendix  E,  the  use  of  the  Acquisition  Range  Method  involves 
a  singular  (i.e. ,  at  least  one  zero  eigenvalue)  covariance  matrix  at  acquisition. 

At  every  sample  time  in  the  baseline  simulation  it  was  noted  that  one  of  the  eigen¬ 
values  of  the  covariance  matrix  was  close  to  zero  in  magnitude  and  not  always 
positive.  This  situation,  however,  had  no  apparent  effect  on  the  performance  of 
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Table  8-1.  Baseline  Configuration  for  One-Way  Doppler  Study 


Key  Parameters 


Measurement  Model 

Acquisition  Range  Reference  (see  Appendix  E) 

Number  of  Stations 

21 

Gravity  Model 

Real  World:  APL  "8-8" 

Filter  World:  h1  and 

Atmosphere  Model 

Real  World:  U.  S.  Standard  1962 

Filter  World:  exponential 

Data  Rate 

1  sample/ 10  seconds 

Measurement  Noise 

50  ft  (l«r) 

Forcing  Noise 

Real  World:  £[  u(t)u(T)]  =  .  01  6(t  -  t)  ft^/sec^ 

(only  on  &y ) 

Filter  World:  none 

Initial  Conditions 


States 

Real  World 

Filter  World 

Standard  Deviation 

xl»x2*x3 

Reference  Orbit 

6700  ft  offset 
(estimated) 

10,000  ft 

x4,x5,x6 

Reference  Orbit 

6.  7  ft  offset 
(estimated) 

10  ft/ sec 

*7 

10,000  ft 

not  estimated 

— 

X8 

5  ft/ sec 

0  ft/ sec 
(estimated) 

5  ft/ sec 

X15 

computed  internally 

computed  internally 

computed  internally 

X16 

.  01  ft/ sec^ 

not  estimated 

— 
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the  filter  on  a  computer  with  a  60-bit  word.  On  the  other  hand,  if  a  much 
smaller  word  length  machine  were  used,  roundoff  could  become  a  significant 
problem. 

8.  3  EFFECT  OF  DATA  RATE 

8.  3.  1  The  baseline  system  with  104  ground  stations  was  examined  with  data 
rates  of  10,  20,  and  30  sec.  The  full  deployment  of  104  ground  stations  was 
chosen  to  minimise  the  effect  of  stations  whose  visibility  period  is  less  than 
30  sec.  The  results  are  summarised  in  Table  8-2. 


Table  8-2.  Effect  of  Data  Rate  on  One-Way  Doppler  RSS  Position  Errors 


Time  from 
Epoch 
(sec) 

10-sec 

Data  Rate 

20- sec 

Data  Rate 

u 

V 

■ 

1 

o 

fO 

Data  Pate 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Star.aard 

Deviation 

(ft) 

0 

11604 

17321 

11604 

17321 

11604 

17321 

1200 

459 

828 

1108 

1114 

1056 

1280 

2400 

145 

153 

479 

218 

264 

261 

3600 

446 

116 

443 

174 

560 

221 

4200 

237 

54 

165 

76 

266 

95 

8.3.2  For  all  rates  examined,  modeling  error  appears  to  become  significant 
after  about  one-half  an  orbit.  The  effect  is,  of  course,  less  prominent  at 
the  30-sec  data  rate.  In  fact,  with  21  ground  stations  and  a  30-sec  data  rate, 
modeling  error  does  not  become  significant  until  the  second  revolution. 

8.3.3  Although  the  overall  effect  of  increased  data  rates  seems  negligible  as 
far  as  accuracy  is  concerned,  a  10-sec  rate  is  preferred  over  the  slower 
rates  since  it  may  permit  use  of  information  from  a  station  whose  visibility  is 
relatively  short  for  a  particular  pass. 
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8.  3.  4  The  10-sec  rate  would  also  tend  to  accelerate  convergence  of  the 
state  for  the  frequency  offset  (xg). 

8.  4  OSCILLATOR  ERRORS 

8.4.  1  For  the  values  considered,  the  effects  of  frequency  offset  drift,  Xj^,  and 
the  white  noise  forcing  on  xy  can  be  igaored  by  the  filter. 

8.  4.  2  In  modeling  the  predicted  range  difference,  the  value  of  Xy  essentially 
cancels  out;  therefore,  the  quantity  xy  is  unobservable.  The  presence  of 
white  noise  on  Xy,  of  course,  leads  to  a  random  walk  on  the  measurement. 
However,  since  the  duration  of  a  pass  is  relatively  short,  the  effect  is  not 
significant  for  the  parameter  values  considered. 

8.  4.  3  The  unintentional  frequency  offset  of  the  oscillators,  Xg,  is  generally 
observable.  For  long-duration  passes  over  a  ground  station,  as  shown  in 
Table  8-3,  almost  all  of  this  offset  can  be  recovered  by  the  end  of  the  pass. 
While  the  offset  rate  may  not  be  fully  estimated  for  short  passes,  its  net 
effect  will,  of  course,  be  small  in  such  cases. 


Table  8-3.  Observability  of  Frequency  Offset  for  One-Way 
Doppler  Measurement 


Time  from  Station 
Acquisition 
(seconds) 

Actual  Frequency 
Offset 
(ft/sec) 

Error  in  Estimated 
Frequency  Offset 
(ft/ sec) 

Standard 

Deviation 

(ft/sec) 

0 

15 

15 

15 

20 

15 

4.2 

9.9 

40 

15 

8.3 

8.9 

60 

15 

10.  1 

8.2 

80 

15 

1.  5 

3.  1 

100 

15 

1.2 

1.7 

120 

15 

1.  2 

1.3  > 

140 

15 

0.9 

0.9 

160 

15 

0.8 

0.  5 

180 

15 

0.2 

0.3 

200 

15 

0.3 

0.2 

220 

15 

0.  1 

0.2 
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8.  4.  4  The  frequency  offset  of  the  oscillators,  Xg,  must  be  estimated  in  the 
filter.  This  is  obvious  from  Table  8-4  where  the  RSS  position  errors  obtained 
when  estimating  or  not  estimating  the  drift  are  summarised. 


Table  8-4.  Effect  of  Frequency  Offset  Estimation  on 
RSS  Position  Error 


Time  of  Acquisition 
and  Termination  of 
Ground  Station 
(sec) 

Ground 

Station 

Number 

5  ft/ sec  Offset 

15  ft/ sec  Offset  ! 

Estimated 

(ft) 

Nbt 

Estimated 

(ft) 

Estimated 

(ft) 

Not 

Estimated 

(ft) 

0 

11600 

11600 

11600 

11600 

940 

19 

19555 

19555 

19555 

19555 

980 

3976 

3696 

4218 

4117 

1180 

4 

4588 

4362 

5046 

4968 

1310 

551 

6760 

3307 

20600 

1560 

3 

1268 

11282 

4998 

34256 

1780 

326 

4650 

887 

11688 

3200 

10 

1285 

23177 

5096 

58779 

3380 

139 

1297 

917 

3991 

8.5  INITIAL  CONVERGENCE  ANALYSIS 

As  shown  in  Tables  8-5  and  8-6,  the  filter  has  difficulty  converging  for  initial 
position  offsets  over  three  n  mi.  These  results  are  similar  to  those  obtained 
with  range- rate  data  in  Sect.  7  and  suggest  the  leed  for  horizon  sensors. 
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Table  8-5,  RSS  Position  Errors  for  Three  n  mi  Initial  Offset 


Time  from 
Epoch 
(sec) 

104  Ground  Stations 

21  Ground  Stations 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

0 

31177 

54772 

31177 

5477  2 

1000 

3078 

403 

19299 

8589 

1300 

1265 

170 

32414 

6807 

1700 

38 

58 

21622 

846 

3300 

175 

29 

12501 

535 

Table  8-6.  RSS  Position  Errors  for  Five  n  mi  Initial  Offset 


Time  from 
Epoch 
(sec) 

104  Ground  Station 

21  Ground  Stations 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

0 

103923 

173205 

103923 

173205 

1000 

763111 

663 

181828 

13533 

1300 

147091 

115 

415801 

1305 

1700 

83546 

67 

751210 

212 

3300 

42811 

28 

3582780 

45 

8.6  STATION  LOCATION  ERRORS 

The  effect  of  errors  in  the  location  of  the  ground  stations  was  investigated. 

In  Tables  8-7  and  8-8  the  results  of  estimating  these  errors  or  ignoring  them 
completely  are  summarized.  It  is  concluded  that  they  should  be  estimated  if 
their  expected  value  is  more  than  100  ft. 


8-11 


Table  8-7.  Effect  of  Station  Location  Errors  at  2400  Sec  on  RSS 
Position  Errors 


Station 

Location 

Errors 

(ft) 

Est 

mated 

Not  Estimated 

RSS 

Error 

(ft) 

RSS  Standard 
Deviation 
(ft) 

RSS 

Error 

(ft) 

RSS  Standard 
Deviation 
(ft) 

0 

264 

261 

264 

261 

50 

251 

271 

269 

261 

200 

354 

543 

456 

261 

400 

693 

698 

831 

261 

Table  8-8.  Effect  of  Station  Location  Errors  at  4200  Sec  on  RSS 
Position  Errors 


Station 

Location 

Errors 

(ft) 

Est 

imated 

Not  Estimated 

RSS 

Error 

(ft) 

RSS  Standard 
Deviation 
(ft) 

RSS 

Error 

(ft) 

RSS  Standard 
Deviation 
(ft) 

0 

265 

95 

265 

95 

50 

323 

106 

344 

95 

200 

396 

286 

583 

95 

400 

601 

412 

903 

95 

8.7  MEASUREMENT  MODELING  APPROACH 

EVALUATION  (SEE  APPENDIX  E) 

8.  7.  1  Three  methods  for  modeling  the  data  were  discussed  in  Section  8.  1.  A 
comparison  of  the  Acquisition  Reference  Range  (ARR)  and  the  Previous  Reference 
Range  (PRR)  methods  for  baseline  conditions  is  summarised  in  Table  8-9.  The 
ARR  method  gives  somewhat  better  navigation  accuracies.  This  result  can  be 
misleading  since  it  is  quite  dependent  on  the  validity  of  the  real-world  data  model 
used.  In  this  model,  it  was  implicitly  assumed  that  measurement  noise  occurs 
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only  upon  readout  of  the  accumulated  range  difference.  From  this,  it  follows 
that  the  assumption  of  white  noise  in  filter  is  strictly  correct  only  for  the  ARR 
method.  Furthermore,  in  the  ARR  method,  the  effect  of  this  readout  noise  is 
clearly  diminished  with  each  successive  measurement.  Thus,  considering  that 
the  assumed  real-world  model  is  quite  unfavorable  to  the  PRR  method,  that 
approach  appears  to  give  remarkably  good  results. 

8.  7.  2  Indeed,  even  if  the  real-world  model  were  correct,  there  are  several 
distinct  advantages  with  the  PRR  method.  The  dimension  of  the  state  vector 
is  smaller,  covariance  initialisation  is  not  required,  and,  if  the  noise  on  Xy 
is  truly  significant,  the  resulting  random  walk  effect  would  considerably  degrade 
the  ARR  method  but  not  the  PRR  method. 

Table  8-9*  Comparison  of  RSS  Position  Errors  Obtained  by  Acquisition 
Acquisition  Range  Reference  and  Previous  Range 
Reference  Methods 


Time  from 
Epoch 
(sec) 

Acquisition  Reference 
Range  Method 

Previous  Reference 
Range  Method 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

0 

11605 

17321 

11605 

17321 

1200 

1056 

1280 

764 

1490 

2400 

264 

261 

254 

242 

4200 

266 

95 

422 

148 

8.  7.  3  The  use  of  instantaneous  range  rate  computed  at  the  middle  of  a  data  col¬ 
lection  interval  as  an  approximate  ion  to  the  measured  average  range  rate  over 
that  interval  was  tested  by  a  series  of  simulations  in  which  the  collection  interval 
was  varied  from  1  to  10  sec.  For  all  cases,  a  new  measurement  was  taken  every 
10  sec.  The  results,  which  are  summarised  in  Table  8-10,  show  that  the 
approximation  ir  only  valid  for  collection  intervals  less  than  5  sec. 
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8.7.4  In  a  separate  study,  the  average  range  rate  and  the  instantaneous 
value  at  the  midpoint  of  the  10  second  data  intervals  were  directly  compared  for 
various  orbital  and  ground  station  configurations.  The  study  showed  that 
the  approximation  in  question  is  quite  good  near  the  time  of  satellite  rise 
or  set,  but  is  generally  quite  poor  (up  to  1  n  mi/sec  error)  at  the  time  of 
closest  approach.  Upon  careful  examination,  this  effect  was  observed  in  the 
simulation  runs  which  tested  this  approximation.  That  is,  the  filter  seemed 
to  perform  adequately  near  each  station  acquisition,  but  eventually  became 
''confused".  Eventually,  of  course,  this  confusion  lead  to  filter  divergence. 


Table  8-10.  RSS  Position  Errors  Obtained  by  Use  of  Instantaneous 
Range  Rate  for  Processing  One-Way  Doppler  Data 


Time  from 
Epoch 
(sec) 

Accumulation  Period  for  Range-Rate  Approximation 
(New  Accumulation  every  10  sec  for  all  cases) 

10  seconds 

5  seconds 

1  second 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

RSS 

Error 

(ft) 

Standard 

Deviation 

(ft) 

0 

11605 

17321 

11605 

17321 

11605 

17321 

1000 

14600 

12894 

15206 

12909 

15331 

12878 

1300 

2726 

6520 

742 

6479 

315 

6546 

1800 

31363 

17349 

7694 

1521 

317 

1642 

3400 

2656 

715 

683 

678 

44 

697  1 
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9.  TWO  ANGLES  TO  KNOWN  LANDMARKS 


9.  1  BASELINE  CONFIGURATION  AND  RESULTS 

Although  extremely  high  accuracy  navigation  was  shown  to  be  possible  with 
range  rate  (or  one-way  doppler)  systems;  there  is  an  admitted  loss  of  auton¬ 
omy.  It  is  possible  that  such  accuracy  would  be  obtained  with  somewhat 
more  autonomy  by  means  of  angular  measurements  to  known  locations.  Such 
measurements  could  be  obtained  in  practice  by  either  optical  tracking  of 
known  surface  features  or  antenna  tracking  of  electromagnetically  radiating 
cooperative  or  uncooperative  sources.  In  this  section  the  generic  problem 
of  navigating  with  such  measurements  is  considered.  The  specific  baseline 
configuration  used  for  this  study  is  described  in  Table  9-1. 

The  results  of  simulating  the  baseline  case  are  shown  graphically  in  Fig.  9-1 
and  summarized  in  Table  9-2.  As  might  be  anticipated,  the  general  behavior 
of  angular  measurements  is  quite  similar  to  range  rate  measurements. 

Indeed,  the  results  seem  somewhat  more  stable  during  the  first  orbit. 

9.2  EFFECT  OF  INCREASED  SENSOR  ACCURACY 

As  indicated  by  comparison  of  Tables  9-2  and  9-3  improving  the  sensor 
accuracy  to  30  arc  sec  has  little  effect  during  the  first  three  landmarks 
while  performance  is  enhanced  by  a  factor  of  two  after  the  third  landmark. 

9.3  EFFECT  OF  HORIZON  SENSOR  AUGMENTATION 

The  results  obtained  by  adding  the  basic  horizon  sensor  configuration  of 
Sec.  5  to  the  30-arc  sec  dual-angle  measurement  system  are  shown  in 
Fig.  9-2  and  summarized  in  Table  9-4.  The  horizon  sensor  system  effects 
a  marked  improvement  in  performance  with  the  first  two  landmarks,  but 
does  not  significantly  change  the  steady-state  errors. 

9.4  EFFECT  OF  INCREASING  LANDMARKS 

The  30-arc  sec  dual-angle  system  augmented  with  horizon  sensors  was  studied 
with  the  full  landmark  deployment  of  Appendix  F.  The  improved  performance 
of  the  system  is  clearly  seen  in  Table  9-5. 
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Table  9-1.  Baseline  Configuration  for  Two  Angles  to  Known 
Landmarks 


Key  Parameters 


Measurement  Model 

Right  ascension  and  declination  to  landmark 
using  an  inertial  coordinate  system  centered 
at  the  navigator 

Number  of  Stations 

21 

Gravity  Model 

Real  world:  APL  "8-8" 

Filter  world:  p  and  J 2 

Atmosphere  Model 

Real  world:  U.S.  Standard  1962 

Filter  world:  Exponential 

Data  Rate 

1  sample/10  seconds 

Measurement  Noise 

60  arc  sec  on  each  angle 

Forcing  Noise 

None 

Minimum  elevation 
angle  for  landmark 
visibility 

5  deg 

Initial  Conditions 


State 

Real  World 

Filter  World 
(all  states  estimated) 

Standard 

Deviation 

Xj  x2  x„  (position) 

Reference  Orbit 

6700  in  each 
direction 

x4  x5  x6  (velocity) 

Reference  Orbit 

6.7  ft/sec  in  each 
direction 

10  ft/sec 
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Figure  9-1.  RSS  Position  Errors  for  Baseline  Dual- Angle  System 


Table  9-2.  RSS  Error  for  Dual-Angle  Baseline  System 


Time  of  Acquisition 
and  Termination  of 
Landmark  (sec) 

Landmark 

Number 

RSS 

Error 

(ft) 

RSS 

Standard 
Deviation  (ft) 

0 

-- 

11604 

17320 

930 

6045 

28938 

990 

10 

2416 

4663 

1170 

4490 

5669 

1320 

4 

1618 

1131 

1530 

2901 

2453 

1800 

3 

309 

286 

3180 

2486 

2247 

3390 

10 

293 

637 

5400 

-- 

1394 

2046 

Table  9-3.  RSS  Errors  for  Dual-Angle  System  Using  30  Arc  Sec 
Angle  Errors 


Time  of  Acquisition 
and  Termination  of 
Landmark  (sec) 

Landmark 

Number 

RSS 

Error 

(ft) 

RSS 

Standard 
Deviation  (ft) 

0 

-- 

11604 

17320 

930 

6044 

28938 

990 

19 

914 

2778 

1170 

3230 

4278 

1320 

4 

661 

623 

1560 

877 

1746 

1800 

3 

112 

144 

3180 

877 

1139  j 

3390 

10 

233 

319 

5400 

-- 

1332 

1028 
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Figure  9-2.  RSS  Position  Errors  for  Dual-Angle  Horizon  Sensor 
Configuration  (Using  30  arcsec  Angle  Errors) 


Table  9-4.  RSS  Position  Errors  for  Dual  Angle-Horizon  Sensor 
System  Using  30  Arc  Sec  Angular  Errors 


Time  of  Landmark 
Acquisition  and 
Termination  (sec) 

Landmark 

Number 

RSS  Position 
Error 
(ft) 

RSS  Standard 
Deviation 
(ft) 

0 

-- 

11604 

17320 

930 

19 

11973 

6963 

990 

1786 

2342 

1170 

1288 

2787 

1320 

4 

797 

563 

1530 

% 

2078 

1193 

1800 

88.4 

142 

3180 

10 

622 

1093 

3390 

196 

318 

5400 

m  m 

805 

1004 

9.5  EFFECT  OF  MULTIPLE  SIMULTANEOUS  LANDMARK 

TRACKING  ON  CONVERGENCE 

From  previous  results*  it  is  clear  that  the  use  of  horizon  sensors  is  a  good 
approach  for  ensuring  convergence  from  large  initial  position  errors.  As  a 
possible  alternative*  the  use  of  simultaneous  tracking  of  more  than  one  land¬ 
mark  is  of  interest.  To  examine  this  approach*  a  configuration  with  three 
different  sets  of  30-arc  sec  angle  trackers  was  devised  and  simulations  were 
performed  for  initial  errors  of  1,  5,  10,  and  20  n  mi.  Perfectly  matched 
dynamic  models  were  used  in  order  to  isolate  the  cause  of  possible  divergence. 

The  results*  as  summarized  in  Table  9-6*  show  that  convergence  is  not 
possible  for  large  initial  errors  if  only  a  filter  algorithm  is  used.  However, 
some  special  initialization  procedure  such  as  triangulation  using  simultaneous 
measurements  to  two  landmarks  could  possibly  be  used  effectively. 
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Table  9-5.  RSS  Position  Error*  for  Dual-Angle  30  Arc  Second  Error- 
Horizon  Sensor  System  with  104  Landmarks 


Time  of  Acquisition 

RSS 

RSS 

1  and  Termination  of 

Number 

Error 

Standard 

Landmarks  (sec) 

(ft) 

Deviation  (ft) 

0 

80 

11604 

17320 

190 

3534 

1789 

200 

72 

649 

469 

390 

433 

448 

870 

63 

1295 

1111 

1150 

45.3 

96.7 

1160 

26 

39.6 

97.3 

1250 

31.3 

102 

1260 

A 

30.9 

104 

1310 

47.9 

104 

1320 

23 

48.4 

106 

1480 

52.2 

73.1 

1490 

25 

48 

73.7 

1560 

74.7 

78.4 

1560 

* 

74.7 

78.4 

1780 

O 

85.1 

67.4 

1790 

24 

89.1 

67.9 

*  1830 

102 

70.9 

2070 

68 

156 

94.1 

2310 

.  109 

77.5 

2310 

70 

110 

77 

2520 

116 

75 

2550 

40 

121 

77.9 

2820 

135 

76.3 

3000 

88 

169 

89 

3300 

130 

78.6 

3310 

10 

129 

78.8 

3390 

132 

81 

3540 

93 

135 

91.9 

3930 

150 

76.6 

3940 

56 

162 

76.8 

4020 

161 

79.8 

4020 

57 

166 

79. 1 

4290 

97.  3 

71.6 

5400 

-- 

203 

138 
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Table  9-6.  RSS  Position  Errors  for  Three  Simultaneous  Dual-Angle  Trackers 


RSS 

Standard 

Deviation 

(ft) 

246410 

731 

1503 

107 

151 

139 

145 

141 

RSS 

Error 

(ft) 

232095 

10510 

10466 

666 

1049 

780 

795 

750 

— 

RSS 

Standard 

Deviation 

(ft) 

173205 

722 

1491 

106 

150 

139 

145 

140 

RSS 

Error 

(ft) 

116047 

3425 

2849 

191 

292 

250 

252 

228 

RSS 

Standard 

Deviation 

(ft) 

86602 

717 

i486 

106 

150 

138 

145 

140 

RSS 
Error 
(ft)  • 

58023 

1438 

1020 

61.8 

98 

97 

96.5 

82.2 

RSS 

Standard 

Deviation 

(ft) 

17320 

705 

1472 

106 

150 

138 

144 

140.6 

RSS 

Error 

(ft) 

11604 

655 

1058 

24.9 

63.5 

42.8 

40 

36.2 

Time 

From 

Epoch 

(sec) 

0 

j  370 

870 

1810 

2090 

2500 

2550 

2700 

10.  OTHER  SYSTEMS 


Since  the  time  available  for  the  current  study  was  quite  limited,  only  those 
systems  which  a  priori  appeared  to  provide  an  optimum  combination  of  accu¬ 
racy,  autonomy,  and  cost  could  be  studied  in  detail.  In  this  section  some  of 
the  systems  which  were  not  simulated  are  discussed  qualitatively. 

10. 1  NAVIGATION  SATELLITES 

For  several  years,  the  concept  of  navigation  by  observing  specially  designed 
satellites  has  received  considerable  attention  (Ref.  44).  The  following  key 
arguments  are  noted: 

10.  1. 1  Since  navigation  satellites  would  essentially  be  known  points  in  inertial 
space,  the  accuracy  performance  should  be  at  least  comparable  to  range  or 
range-rate  observations  to  ground  stations.  Performance  with  angle  measure¬ 
ments  would  depend  strongly  on  the  relative  distance  between  the  user  and  the 
navigation  satellite. 

10. 1.2  A  very  attractive  property  of  satellites  is  their  wide  range  of  visibility 
as  compared  to  ground  stations.  In  particular, with  synchronous  navigation 
satellites,  the  coverage  area  is  practically  a  hemisphere  for  low-altitude  navi¬ 
gators.  Indeed,  the  coverage  area  grows  even  larger  as  the  navigator's  altitude 
increases.  Global  coverage  could  easily  be  obtained  with  a  relatively  small 
number  of  satellites  if  they  were  properly  located. 

10.  1.3  In  order  to  use  satellites  for  navigational  observations,  their  ephem- 
erides  must  be  maintained.  This  implies  the  utilization  of  a  ground-based 
tracking  network  and  complex.  Clearly,  such  a  system  cannot  inherently  pro¬ 
vide  a  high  degree  of  autonomy  to  a  user. 
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10.2 


EJECTED  PROBE 


The  concept  of  orbit  navigation  by  observing  a  secondary  "mini"  satellite  which 
has  been  ejected  from  the  primary  vehicle  has  recently  been  investigated 
(Ref.  38).  The  relative  motions  of  the  two  satellites  when  coupled  with  the 
dynamical  models  for  both  can  theoretically  provide  enough  information  for 
navigation  in  inertial  space.  However,  the  accuracies  reported  in  the  reference 
are  not  very  encouraging.  Even  if  all  model  errors  were  neglected,  accuracies 
of  only  3633,  1739,  and  1358  feet  after  1,  2,  and  3  orbits,  respectively,  were 
indicated.  Nevertheless,  because  of  the  inherent  high  degree  of  autonomy,  the 
concept  should  be  analyzed  further  before  it  is  completely  eliminated  from  con¬ 
sideration. 

10.3  STAR-HORIZON  MEASUREMENTS 

A  possible  approach  to  autonomous  navigation  could  be  based  on  observing  the 
angle  between  a  known  star  and  the  horizon  of  the  earth  as  measured  in  the 
plane  defined  by  the  navigator,  the  star,  and  the  center  of  the  earth.  The  infor¬ 
mation  contained  in  this  so-called  "star -horizon"  or  "star-elevation"  measure¬ 
ment  is  essentially  the  same  as  that  of  the  horizon  sensor  measurement  des¬ 
cribed  in  Section  5;  that  is,  since  horizon  sensor  measurements  inherently 
require  precise  attitude  information  and  since  such  data  can  only  be  achieved 
by  using  star  observations,  the  horizon  sensors  effectively  measure  angles 
between  the  stars  and  the  horizon.  The  basic  difference  lies  in  the  mechaniza¬ 
tion.  While  stars  must  be  acquired  by  the  star- elevation  sensor,  there  is  no 
need  for  high  precision  attitude  information.  Another  important  distinction 
between  the  two  types  of  measurements  resides  in  the  fact  that  horizon  sensor 
observations  are  always  made  in  the  planes  of  the  sensor  heads,  while  the 
star-horizon  measurement  plane  varies  with  the  stars  being  used. 

It  was  shown  in  Section  5  that  deterministic  navigation  was  theoretically 
possible  using  horizon  sensors.  It  can  be  shown  that  three  properly  oriented 
star -elevation  measurements  could  also  provide  such  a  capability.  Thus  no 
significant  difference  in  performance  should  be  expected  between  these  two 
measurement  types. 
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10.4 


STAR  OCCULATION 


Another  approach  to  autonomous  navigation  involves  measurement  of  the  time 
at  which  selected  stars  pass  behind  the  earth's  horizon.  While  this  star  occu- 
lation  is  physically  quite  different  from  the  star -elevation  observation  des¬ 
cribed  in  Section  10.  3,  the  navigation  information  contained  in  both  measure¬ 
ments  is  basically  the  same.  It  can  be  shown  (Ref.  3)  that  both  measurements 
provide  position  information  normal  to  the  line  of  sight  to  the  horizon  in  a  plane 
defined  by  the  navigator,  the  star,  and  the  center  of  the  earth. 

The  major  difference  between  these  two  sensors  lies  in  the  phenomena 
measured;  that  is,  horizon  sensors  depend  on  the  carbon  dioxide  absorption 
layer  of  the  atmosphere,  while  the  occulation  measurement  would  depend  on  the 
instantaneous  atmospheric  refraction  and  absorption  of  light  at  the  time  of  mea¬ 
surement.  Notwithstanding  these  model  differences,  however,  the  convergence 
properties  of  an  occulation  system  should  be  similar  to  the  horizon  sensing 
system  already  discussed.  Specification  of  the  accuracy  attainable  would,  of 
course,  depend  on  detailed  model  analysis  and  simulation. 
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1  i .  COMPUTER  WORDLENGTH  EFFECTS 


Numerical  roundoff  difficulties  in  recursive  filtering  have  been  widely  studied 
and  publicized  in  recent  years.  Although  many  techniques  for  alleviating 
these  difficulties  have  been  suggested,  all  are  ultimately  limited  by  the  basic 
wordlength  of  the  machine  used  to  process  data. 

In  order  to  determine  the  sensitivity  of  satellite  navigation  to  such  computer 
characteristics,  two  baseline  configurations  with  supposedly  different  numerical 
properties  were  simulated  with  a  variety  of  wordlengths.  Square  root  filtering 
was  used  in  all  cases  to  enhance  the  overall  performance.  The  unknown  land¬ 
mark  horizon  sensor  baseline  configuration  was  used  as  one  standard  because 
of  the  relatively  high  dimensional  state  vector  (14)  involved,  and  also 
because  unknown  landmark  tracking  has  been  reported  t~  be  particularly 
sensitive  to  numerical  errors.  To  maintain  perspective,  the  straight¬ 
forward  case  of  range  rate  to  known  landmark  was  chosen  as  the  other 
basic  example.  This  configuration  required  only  six  states  for  navigation 
and  no  peculiar  numerical  problems  were  noted. 

All  simulations  were  performed  on  a  Control  Data  6600  60-bit  wordlength 
machine  by  means  of  special  software  which  could  truncate  the  results  of 
selected  arithmetic  operations.  To  compensate  for  optimum  coding  that  would 
undoubtedly  be  provided  for  operational  navigation  software,  only  those  calcu¬ 
lations  directly  associated  with  matrix  multiplication,  addition,  or  subtraction 
were  actually  subjected  to  truncation.  That  is,  matrix  inversion,  square- root 
extraction,  integration,  and  partial  derivative  computations  were  all  performed 
with  the  full  60  bit  word  of  the  6600  computer.  The  results  reported  should 
be  highly  optimistic,  since  failure  of  a  simulated  truncation  under  these 
conditions  would  surely  imply  failure  if  all  operations  were  truncated. 

Under  the  above  assumptions,  the  baseline  configurations  chosen  were  simulated 
with  60,  30,  20  and  15  bit  words.  The  results  obtained  are  summarized  in 
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Table*  11-1  and  11-2.  Comparison  of  these  tables  shows  that  in  spite  of  the 
advertised  numerical  differences,  the  configurations  behaved  remarkably 
similar  with  respect  to  wordlength  limitation.  It  appears  that  a  30-bit 
machine  would  be  adequate  for  both  systems,  while  use  of  only  20  bits  could 
result  in  considerable  numerical  difficulty.  Since  the  study  was  admittedly 
highly  optimistic,  it  is  not  likely  that  machines  with  less  than  20  bits  could 
be  effective  for  navigation. 

The  authors  would  like  to  acknowledge  the  assistance  of  Dr.  H.  J.  Wertz, 
who  prepared  the  special  software  that  was  used  to  simulate  the  various 
computer  wordlengths. 


Table  11-1.  Effect  of  Computer  Wordlength  With  Baseline  Unknown  Landmark 
Tracking  -  Horizon  Sensor  Configuration 
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Table  11-2.  Effect  of  Computer  Word  Length  With  Baseline  Range  Rate 
to  Known  Landmarks  Configuration 
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Standard  '■ 
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(ft) 
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Error 

(ft) 

10407 

26035 

69823 

155941 

215481 

RSS 

Standard 

Deviation 

(ft) 

17320 

619 

118 

104 

68 

RSS 

Error 

(ft) 

11574 

1244 

2051 

4994 

6558 

RSS 

Standard 

Deviation 

(ft) 

17320 

619 

118 

104 

68 

RSS 

Error 

(ft) 

11604 

573 

122 

445 

347 

RSS 

Standard 

De  viation 
(ft) 

17320 

619 

118 

104 

68 

RSS 

Error 

(ft) 

11604 

573 

121 

441 

341 

Time  from 
Epoch  (sec) 

0 

1000 

2090 

3010 

4000 

12.  PRIMARY  CONCLUSIONS 


12. 1  Horizon  senior*  appear  to  be  a  promising  candidate  for  inclusion  in 
hybrid  autonomous  navigation  systems  in  order  to  provide  recovery  from 
large  off-nominal  conditions,  bound  maximum  errors,  and  accelerate  filter 
convergence.  On  the  other  hand,  the  use  of  horizon  sensors  alone  can  only 
provide  coarse  navigation  accuracy. 

12.2  The  convergence  properties  of  horizon  sensor  systems  are  enhanced  by 
inclusion  of  the  measurement  of  the  angle  subtended  by  the  earth's  disk. 

12.  3  The  ability  to  recover  from  large  initial  errors  is  not  a  natural  charac¬ 
teristic  of  many  navigation  configurations.  The  associated  problem  generally 
requires  special  attention. 

12.4  Unknown  landmark  tracking,  while  attractive  from  an  autonomy  point  of 
view,  provides  only  moderate  accuracy  with  oscillatory  behavior  unless  a 
high  accuracy  sensor  and  a  relatively  uniform  distribution  of  landmarks  with 
small  altitude  uncertainties  are  ensured.  To  achieve  this  uniform  distribution 
of  landmarks,  the  night-day  and  cloud-cover  problems  must  be  solved. 

12.5  Systems  utilizing  optical  or  electronic  measurements  (e.g.,  range, 
Doppler,  or  angular  measurements)  to  known  landmarks  or  satellites  can  be 
used  to  provide  high  navigation  accuracy.  The  position  accuracy  of  the  land¬ 
marks  or  satellites  is  generally  more  critical  than  the  choice  of  the  sensor 
type  used. 

12.6  Even  in  the  presence  of  time-varying  oscillator  model  errors,  one-way 
Doppler  systems  can  yield  excellent  navigation  accuracy.  Because  of  the 
short  duration  of  a  station  pass,  the  only  oscillator  error  that  need  be  modeled 
in  the  filter  is  a  constant  frequency  offset  error  between  the  ground  and  user 
oscillators.  This  error  is  generally  observable  for  those  stations  close  to  the 
orbit  plane. 
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12.7  Computer*  for  satellite  navigation  with  30-bit  word  lengths  appear  to 
be  adequate  for  avoiding  filter  divergence  due  to  numerical  round-off.  Based 
on  the  highly  optimistic  study  performed,  it  is  not  likely  that  machines  with 
less  than  20  bit  words  could  be  effective  for  navigation. 


APPENDIX  A 


REAL  WORLD  DYNAMIC  MODELS 

A.  I  GEOPOTENTIAL 


In  general,  geopotential  forces  can  be  found  by  computing  gradients  of 
truncated  versions  of  the  function 
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+  V  Pm(sin  *)(C  cos  m\  +  S  sin  mK) 
/  >  n  nm  nm 
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where 

(i  is  the  product  of  the  universal  gravitational  constant  and  the  mass 
of  the  earth 

r,  p,  and  X  are  the  vehicle  distance  from  the  center  of  the  earth,  the 
vehicle  geocentric  latitude,  and  the  vehicle  longitude 

Re  is  the  earth's  equatorial  radius 

PQ  are  the  Legendre  polynomials  of  the  first  kind  of  degree  n 

P™  are  associated  Legendre  functions  of  the  first  kind  of  degree  n 
and  order  m 

C  are  zonal  harmonics 
n 

C  S  _  are  nonzonal  (sectoral  and  tesseral)  harmonics 
nm  nm 

The  specific  truncation  and  appropriate  coefficients  for  the  model  used  in 
SNAP  were  adopted  from  John  Hopkins  Applied  Physics  Laboratories  (APL) 
"8-8"  model  (Refs.  17  and  30).  The  values  used  are  given  in  Table  A-l 
and  are  based  on  the  following  normalization  formulas. 


A-l 


ATMOSPHERE 


nm 


=  KC 


nm 


3  =  KS 

nm  nm 


A.  2 


The  atmosphere  model  used  was  the  United  States  Standard  Atmosphere, 
1962  (Ref.  43).  A  ballistic  coefficient  of  0.  02  ft^/lb  was  used. 
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Table  A-l.  Values  for  C  and  S  for  APL  "8-8"  Geopotential  Model 
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APPENDIX  B 


FILTER  WORLD  DYNAMIC  MODELS 


B.  1  GEOPOTENTIAL 

Filter  computations  were  based  on  the  simple  geopotential  function 


U  =  [p/r  -  (pR^/r3)]j2P2(sin  *) 


where  the  notation  is  identical  to  that  of  Appendix  A  and  J2  ■  -  C2< 

B.2  ATMOSPHERE 

The  densities  used  were  computed  from  a  formula  developed  by  Schusterman 
(Ref.  40);  that  is.,  the  density  p  at  an  altitude  z  above  the  earth  was  found 
according  to 


**  ~  ^0  [a  +  b(z 


where 

a  =  16.  4  km 
b  =0.232 

Pq  =  3.951  x  IQ’12  gm /cm3 
Zq  =  140  km 

A  ballistic  coefficient  of  0.  02  ft2 /lb  was  used. 
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APPENDIX  C 


REFERENCE  ORBIT 


Four  representations  for  the  initial  conditions  of  the  reference  orbit  are 
given  below.  Key  fundamental  constants  are: 

Epoch  date  =  2230  GMT,  29  Feb.  1980 

p  =  1.4076538841X1016  ft3/sec2 

Earth  radius  =  20,  925,  738  ft 

Earth  rotation  rate  =  4.  1780742X10” 3  deg/ sec 

1.  Earth- Centered  Inertial  (ECI).  Based  on  an  earth-centered  coor- 
dinate  system  with  an  equatorial  x-y  plane  whose  x  and  z  axes  are 
directed  to  the  mean  equinox  and  north  pole,  respectively. 

x  =  -1.  59387494563X107  ft 
y  =  1.  51831714534X107  ft 

z=  3.  84199431557X10" 15  ft 
x  =  6.  1644 1282935X1 03  ft/sec 
y  =  6.  00674668844X1 03  ft/sec 
i  =  2.  36312409387X104  ft/sec 

2.  Spherical 

Right  ascension  (positive  east  from  x-axis)  = 

136.  390753744  deg 

Geocentric  declination  =  1.  0X10"2®  deg 

Angle  between  velocity  vector  and  geocentric 
vertical  =  90.7298014230  deg 

Azimuth  of  velocity  vector  from  north  (Az)  =  340  deg 

Magnitude  of  position  vector  (R)  =  2.  20130059195X10’7  ft 

Magnitude  of  velocity  vector  (V)  =  2.  5149881 5029X104  ft/sec 


3. 


Classical 


Semi-major  axis  =  2.  1776394691 6X1 07  ft 
Eccentricity  =  1.67413647743X10“^ 

Inclination  =  110  deg 

Right  ascension  of  ascending  node  =  136.  3907S3744  deg 

Argument  of  perigee  =  131.  193639651  deg 

Time  of  last  perigee  passage  =  -3442.  155976  sec  from  epoch 

4.  Geographic 

Geodetic  latitude  =  0  deg 

Geodetic  longitude  =  360  deg 

Angle  between  velocity  vector  and  geocentric 
horizontal  =  -0.  7298014230  deg 

Az  =  340  deg 

R  =  2.2013005915X107  ft 

V  =  2.  51498815029X104  ft/sec 
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HORIZON  ERROR  MODEL 


As  described  in  Sect.  5.1,  the  measurements  considered  to  be  made  by 
horizon  sensors  are  the  direction  of  the  local  vertical  defined  by  0j  and  0j 
and  the  angle  subtended  by  the  widest  diameter  of  the  earth's  apparent  disk 
02*  Actually,  0J,  0j,  and  cannot  be  measured  directly  but  are  inferred 
from  the  directions  to  the  horizons.  The  plane  defined  by  the  z^-axis  (see  Fig¬ 
ures  5-1,  5-2,  and  the  apparent  center  of  the  earth  is  shown  in  Figure  D-l.  The 
angles  and  are  the  expected  measurements  to  the  nominal  CO2  absorp¬ 
tion  layer,  while  the  angles  a-j  and  are  the  instantaneous  measurements  to 
the  true  CO2  absorption  layer.  The  measurements  can  be  thus  expressed  as 


“1  ■  «i„  -  6«) +  »i 

*2  *  "2n  +  ‘"2  +  u2 


(D-l) 


where  6orj  and  602  *re  the  errors  due  to  the  uncertainty  in  the  altitude  of  CO2 
absorption  layer  at  the  point  of  tangency,  while  Uj  and  are  white  noise 
processes  representing  sensor  noise. 

The  horizon  errors  have  been  statistically  described  by  autocorrelation 
functions.  Simulation  of  this  type  of  error  is  readily  achieved  by  means  of 
first-order  Markov  processes.  Before  this  can  be  done,  however,  the  com¬ 
plete  autocovariance  matrix  must  be  known.  Me  Arthur  (Ref.  34)  indicates 
that  the  errors  in  the  altitude  of  the  CO2  absorption  layer  are  exponentially 

correlated  in  space  and  time  with  a  great  circle  space  constant  of  2500  n  mi 

2  2 

and  a  time  constant  of  10  days.  The  variance  is  reported  to  be  2887  ft  . 
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If  an  altitude  of  100  n  mi  is  considered,  the  impact  of  this  error  on  the 
measurement  can  be  computed  by  using  the  distance  to  the  horizon  (916.2  n  mi) 
to  compute  the  equivalent  angular  uncertainty 


*«  <D-2> 

Further,  since  the  vehicle’s  ground  trace  travels  along  the  surface  of  the 
earth  at  an  average  rate  of  4  n  mi /sec,  an  equivalent  time  constant  of  the 
space  correlation  can  be  computed  as: 

t  =  2500  n  mi  (1  sec/4  n  mi)  =  625  sec  (D-3) 

c 

An  approximation  for  the  autocovariances  of  60j  and  6a 2  can  now  he  written 
as 


[Sa^tJSa^t  +  x)]  =  <Tj  exp[  -  (t  -  T)/tc] 
[6a2(t)6a2(t  +  T)]  =  exp[  -  (t  -  T)/tc] 


(D-4) 


where  -  0.02971^  deg^  and  tc  =  625  sec.  Note  that  the  10  day  time 

constant  has  been  neglected  for  this  approximation  since  it  is  much  longer 
than  the  period  of  an  orbit,  which  was  the  duration  of  a  typical  simulated 
flight . 


To  determine  the  cross  correlation  between  6a^  and  6a2>  one  can  consider 
the  great  circle  separation  distance  between  the  points  of  tangency  of  the  two 
lines  to  the  horizons.  For  the  orbit  considered  in  this  study,  the  separation 
dirtance,  AD,  varies  from  1423  to  2547  n  mi.  Evaluating  the  appropriate 
exponentials  results  in  an  average  correlation  coefficient  of  approximately 


p  =  exp(- AD/2500)  «  0.  5 


(D-5) 


Thus 


[  6<»1  (t )6ar2(t  +  T)]  =  po,1o-2  exp[  -  (t  -  r)/tc] 

=  (0.  5)(0. 0297)2  exp[  -  (t  -  r  )/t  ] 

C 


(D-6) 


At  this  point,  enough  information  about  6a  ^  and  is  known  to  simulate  their 
influences  on  the  measurements  by  means  of  Markov  processes.  From  Figure 
D-l,  it  can  readily  be  shown  tfcr.t  0j  and  6^  are  linearity  related  to  a ^  and  a ^ 
by  the  transformation 


where 


(D-7) 


From  Eq.  (D-l)  and  (D-7)  the  autocovariance  of  the  horizon  errors  in 
0- space  can  be  found  to  be 


£[60(t)6eT(t  +  t)]  =  a£[  6a(t)6aT(t  +  t)]AT  (D-8) 
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where 


66  =  A6o 


6  a  = 


66  = 


Expanding  Eq.  (D-8), 


[  66(t)66T  (t  +  t)] 


2  2  + 

e 

-1  lj  -6a2(t)6aj(t  + 


t)  -6«j(t)6a2(t  +  t)1  j"  j  j 


6o2(t)6a2(t  +  t)J  [_-l  i 


1  1 


r  h  j  i  ai  -*iv  i  z 

_  SLjljI 

t  7 

L  c  J  -1  1  -<r_<r1p  <rf  -i  i 


=  exp 


zV 


1  1 
7  2 


=  exp 


n 

.  (t  -  t)  * 

L  J 


+,r2  ’  2P°’i°’z) 


1/2  2 

irz  -  *1 


2  2 
"I+'z  + 


2P<ri°’2J 


Now  since  tr ^  =  <r^  from  Eq.  (D-4),  60  and  662  must  be  independent. 


(D-9) 


t 


( 


Further,  if  the  value  of  p  s  0.5  is  substituted,  Eq.  (D-9)  reduce*  to 


£  [  68(t)69 1  (t  +  t)]  =  exp 


.  (L-  t) 


•i 

V\ 


3<r 


(D-10) 


The  horizon  error*  for  0.  and  6,  can  thus  be  simulated  by  two  independent 

1  *2  2 

processes  whose  variances  are  (Tj/4  and  3<Tj,  respectively,  and  whose 
equivalent  time  constant  is  tc>  By  an  identical  analysis,  the  horizon  errors 
for  0'  can  be  simulated  by  a  process  with  the  same  statistics  as  the  process 
for  0j.  Although  there  may  be  correlation  between  0j  and  0j,  they  were 
assumed  to  be  independent.  This  assumption  was  made  because  sufficient 
data  were  not  available  to  realistically  simulate  these  correlations. 

In  summary  then,  the  contribution  of  the  horizon  errors  can  be  modeled  by 
three  independent  Markov  processes.  From  Eq.  (D-10)  formulas  for  genera¬ 
ting  the  instantaneous  values  of  these  errors  can  readily  be  found.  Specifi¬ 
cally,  if  Xj,  ^2*  .  .  .  ,  Xg  are  reserved  for  the  inertial  position  and  velocity 
states  of  the  navigator,  and  x?,  Xg,  x^  are  used  to  denote  the  states  for  these 
biases,  it  can  be  shown  that 

*?  =  -x?/625  +1^,  £^(0)]=  <r2/4  =  .014852 

Xg  =  -x8/625  +  ug,  ^c2(0)J=  <t2/4  =  .014852  (D-ll) 

*9  =  -x9/625  +u9,  ^c2(0)]=  3<r2  =  .051452 

where  u^,  Ug.  are  zero  mean  white  noise  processes  whose  variances  are 


^[u7(t)u7(r)]=?|?(.OI485)26(t  -  T) 

£(u8(t)ug(T)]  =  (•°i4S5)26(t  -  r)  (D- 12) 

£[u9(t)u9(r)]  =  -gfj  (. 05145)26(t  -  T) 

and  6 (t  -  t)  i»  the  Dirac  delta  function. 
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APPENDIX  E 


ONE-WAY  DOPPLER  MODELING 


A*  Indicated  in  Sect.  8.1*  the  Kalman  filter  formulation  cannot  directly 
accommodate  meaeurementa  which  are  functions  of  the  states  at  different 
times.  The  three  methods  of  obtaining  compatibility  which  were  used  in 
this  study  are  described  in  detail  here. 

E.l  ACQUISITION  RANGE  REFERENCE  METHOD 

(P.  Soule.  Ref.  41) 

In  the  acquisition  range  reference  method  the  time  t^  in  Eq.  (8-4)  is  fixed 
at  tg*  the  time  of  acquisition  of  a  ground  station*  and  the  corresponding 
range,  p(tg),  is  made  a  state  variable  x^.  Then  with  measurements  being 
taken  attj,  t2»  •  .  .*  t^as  the  satellite  passes  over  the  ground  station, 
the  data 


yt  s  p(*g)  +  -  x|5  "  W  +  vi 

y2  *  p(t2)  +  -  x15  -  x^tg)  +  v2 

•  •  •  • 

(E-l) 

•  «  i  • 

*  P<tk)  +  W  *  *15  *  W  +  vk 


readily  conform  to  the  usual  Kalman  filter  formulation.  It  is  noted  that, 
since  the  state  p(tg)  is  perfectly  correlated  with  the  position  states  of  the 
satellite  and  the  ground  station  at  tg,  the  correct  state  covariance  matrix 
at  tg  will  be  positive  semidefinite.  In  theory*  however,  there  should  be 
no  difficulty  processing  the  data  since,  as  Kalman  and  Bucy  (Ref.  29)  point 


E-i 


out,  the  Riccati  equation  will  converge  to  a  unique  solution  even  with  a  positive 
semidefinite  initial  condition  under  very  general  conditions. 


The  correct  covariance  matrix  at  acquisition  is  constructed  from  the  pro¬ 
pagated  state  covariance  matrix  by  first  zeroing  out  all  rows  and  columns 
beyond  six.  The  variances  for  states  Xjt  Xg,  .  .  .  ,  x^,  Xjg  are  then  reset 
to  their  initial  values,  forming  a  temporary  covariance  matrix, P^.  Now,  as 
indicated  above,  the  desired  covariance  matrix  must  contain  terms  to  account 
for  the  correlation  that  exists  between  the  navigator's  position  and  the  range 
at  acquisition,  Xjg.  This  matrix,  p£,can  be  computed  directly  from  the 
relationship 

P+  =  F  P’  FT 

k  k 

where  the  transformation  matrix  F  is  defined  by 

sk  =  rsk  (E-2) 

with 

x^  =  the  state  vector  immediately  after  acquisition 


x^  =  the  state  vector  immediately  before  acquisition 


F  = 


»*k 


14X14 


'14X2 


a* 


15 


as 


15 


8* 


15 


a*. 


ax 


as 


16 


1X14 


(E-3) 


and 

*15  =  p*fc0*  =  [(*1  ‘  *9)  +  (*2  "  *lo)  +  (*3  *  xll)  ] 


i.e.,  F  is  the  16  X  16  identity  matrix  with  the  15th  row  reset  as  indicated. 
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It  should  be  noted  that  this  type  of  re- initialization  technique  can  ou  useful 
for  other  sensors  as  well.  In  particular,  it  can  be  used  when  a  new  state 
variable,  which  is  introduced  at  an  arbitrary  sample  time,  has  a  functional 
dependence  on  other  estimated  state  variables  that  are  continuously  included 
in  the  dynamical  model.  Specifically,  this  same  technique  has  proved  very 
useful  in  the  initialization  of  the  covariance  matrix  for  unknown  landmarks, 
especially  when  the  state  vector  included  several  biases  from  multiple  sensor 
types.  When  compared  to  the  conventional  methods  of  calculating  the 
covariance  terms  in  previous  studies  (Refs.  4,  5,  and  17),  the  technique 
suggested  here  appears  to  be  much  more  systematic,  thereby  preventing 
inadvertent  omissions  of  necessary  terms. 

E.2  PREVIOUS  RANGE  REFERENCE  METHOD 

(I.  A,  Gura) 

In  the  previous  range  reference  method,  the  necessity  for  the  state  xJ5  and 
the  accompanying  covariance  initialization  process  can  be  eliminated  by 
choosing  t^  in  Eq.  (8-3)  to  be  the  time  of  the  previous  measurement.  Thus, 

yt  s  P(tj)  -  p(tQ)  +  x?(tt)  -  x^)  +  Vj  -  vQ 
Yz  ~  P<12*  ’  Pftj)  +  x7<t2)  “  atyCtj)  +  v2  ’  vi 

•  s  • 

(E-5) 

•  •  • 

yk  =  -  P(tk.,)  +  x7(tk)  -  x^.,)  +  vk  -  Vl 

In  order  to  show  that  this  formulation  can  be  accommodated  by  the  Kalman 
filter  algorithm,  note  that  the  predicted  measurement  can  be  expressed  as: 

(E-6) 


i 
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where  the  estimated  state  at  based  on  data  taken  to  tj^j, 

and  i*  the  estimated  state  at  t^  based  on  data  taken  to  t^j.  That  is, 

*k/k-l  a  g**k-i/k-i*  *E“7* 

where  g  (  )  is  the  vector  function  symbolising  a  one  step  integration  of  the 
state  vector. 

Now, 

*k-l/k-l  *  g  1(4k/k-i)  (E*8) 

where  the  notation  g"*  (  )  is  used  to  indicate  a  one  step  backwards  integration. 
Thus  Eq.  (E-6)  becomes 

•  P(*k/k.,)  -  ■>[«'1<*k/k-i>]  <E-’> 

which  is  clearly  within  the  domain  of  the  filter  algorithm.  From  Eq.  (E-9) 
the  measurement  partials  can  readily  be  derived  to  be 

£  =  8P-^~  -  *  8P-^k"|i>  -  fe  +  l.k)  (E-10) 

where  $(k  +  l,k)  is  the  transition  matrix  for  the  state  vector  in  question. 
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Not*  from  Eq.  (E-5)  that  tho  measurement  noise  is  not  white  as  is  required 
for  the  Kalman  filter.  That  is,  if 


*<vk  ‘  vk-l>  (vj  ’  vj-i> 


2R  j«k 

-R  j  =k - 1  or  k+1 
0  for  all  other  j 


(E-  1ft) 


For  the  problem  at  hand  however*  R  is  quite  small  relative  to  other  variances, 
and  this  effect  should  be  Insignificant  if  ignored.  This  contention  was  tested 
in  with  the  simulation  program  and  found  to  be  quite  valid  (See  Section  8.7). 

E.3  RANGE  RATE  AT  MIDPOINT 

Consider  Eq.  (E-5).  A  simple  approach  to  data  processing  would  be  to  model 
the  data  using  the  range  rate  at  the  midpoint  of  the  interval  in  question. 

That  is,  the  data  actually  processed  becomes 

?k *  *k  -  ^k-i  (E-l2> 


while  the  predicted  measurement  in  the  filter  is  the  instantaneous  range 
rate  at  (t^  +  t^  j)/2,  multiplied  by  the  interval  length  t^  -  t^  j.  The 
measurement  equation  thus  becomes 

*  w«k +  *k-i>«] +  xs>  “k  -  W +  vk  •  Vi  <E-‘3> 


Although  such  modeling  may  be  intuitively  appealing,  the  degradation  in 
performance  was  shown  to  be  considerable  when  t^  -  t^_j  was  more  than 
five  seconds  (see  Section  8.  8). 
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LOCATIONS  FOR  KNOWN  LANDMARKS  (GROUND  STATIONS) 


GEOCENTRIC 

EAST 

STATION 

LATITUDE 

LONGITUDE 

1 

-23.  2 

314.  1 

2 

40.7 

141. 3 

3 

61.  2 

210.  1 

4 

76.  5 

291. 2 

5 

-77.8 

166.6 

6 

-4.  6 

55.4 

7 

32.2 

253.2 

8 

51.  1 

358.  9 

9 

39.  1 

283.  1 

10 

-34.6 

138.6 

11 

-25.9 

28.  3 

12 

-14.3 

189.2 

13 

14.9 

120.0 

14 

30.3 

262.  2 

15 

34.5 

69.2 

16 

34.  1 

240.8 

17 

19.2 

166.  6 

18 

18.7 

98.9 

19 

63.9 

337.4 

20 

47.3 

306.0 

21 

35.  1 

33.3 

22 

66.5 

210.0 

23 

69.  5 

260.0 

24 

60.0 

195.0 

25 

60.0 

240.0 
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GEOCENTRIC 

EAST 

STATION 

LATITUDE 

LONGITU 

26 

60.  0 

285.  0 

27 

60.  0 

0.  0 

28 

45.  0 

150.  0 

29 

45.  0 

255.  0 

30 

45.  0 

285.  0 

31 

30.  0 

135.  0 

32 

30.  0 

36.  0 

33 

15.  0 

150.  0 

34 

15.  0 

210.  0 

35 

15.  0 

330.  0 

36 

15.  0 

15.0 

37 

15.  0 

45.  0 

38 

15.  0 

90.  0 

39 

0.  0 

135.  0 

40 

0.  0 

165.  0 

41 

0.  0 

195.  0 

42 

0.  0 

285.  0 

43 

0.  0 

315.  0 

44 

0.  0 

15.  0 

45 

0.  0 

45.  0 

46 

0.  0 

105.  0 

47 

-15.0 

135.  0 

48 

-15.  0 

300.  0 

49 

-15.  0 

15.  0 

50 

-15.  0 

45.  0 

51 

-30.  0 

135.  0 

52 

-30.  0 

180.  0 

53 

-30.  0 

300.  0 

54 

-30.  0 

15.0 

55 

-45.  0 

172.  0 

i 
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GEOCENTRIC 

EAST 

STATION 

LATITUDE 

LONGITUDE 

56 

-66.  5 

105.  0 

57 

-66.  5 

60.  0 

58 

-66.5 

90.  0 

59 

66.  5 

225.  0 

60 

66.  5 

15.0 

61 

60.  0 

225.  0 

62 

60.  0 

255.  0 

63 

60.  0 

315.  0 

64 

60.  0 

15.  0 

65 

45.  0 

240.  0 

66 

45.  0 

270.  0 

67 

45.  0 

0.  0 

68 

30.  0 

180.  0 

69 

30.  0 

75.  0 

70 

15.  0 

180.  0 

71 

15.  0 

300.  0 

72 

15.  0 

0.  0 

73 

15.  0 

30.  0 

74 

15.  0 

75. 0 

75 

0.  0 

120.  0 

76 

0.  0 

150.  0 

77 

0.  0 

180.  0 

78 

0.  0 

285. 0 

79 

0.  0 

300.  0 

80 

0.  0 

0.  0 

81 

0.  0 

30.  0 

82 

0.  0 

75.  0 

83 

-15.  0 

120.  0 

84 

-15.  0 

285.  0 

85 

-15.  0 

315.  0 
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GEOCENTRIC 

EAST 

STATION 

LATITUDE 

LONGITUDE 

86 

-15.0 

30.0 

87 

-30.0 

120.0 

88 

-30.0 

150.0 

89 

-30.0 

285.0 

90 

-30.0 

315.0 

91 

-30.0 

30.0 

92 

-45.0 

285.0 

93 

-66.5 

120.  0 

94 

-66.5 

300.0 

95 

10.0 

250.0 

96 

-15.0 

210.0 

97 

-23.5 

225.0 

98 

-23.  5 

240.0 

99 

-23.  5 

255.0 

100 

-37.0 

210.0 

101 

-47.0 

’7.5 

102 

-47.0 

52.5 

103 

-47.0 

67.5 

104 

-53.  0 

330.0 
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